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Abstract 



A discrete-time quantum walk on a graph is the repeated apphcation of a unitary evolu- 
tion operator to a Hilbert space corresponding to the graph. Hitting times for discrete 
quantum walks on graphs give an average time before the walk reaches an ending con- 
dition. We derive an expression for hitting time using superoperators, and numerically 
evaluate it for the walk on the hypercube for various coins and decoherence models. The 
hitting time for a classical random walk on a connected graph will always be finite. We 
show that, by contrast, quantum walks can have infinite hitting times for some initial 
states. We seek criteria to determine if a given walk on a graph will have infinite hit- 
ting times, and find a sufficient condition, which for discrete time quantum walks is that 
the degeneracy of the evolution operator be greater than the degree of the graph. The 
phenomenon of infinite hitting times is in general a consequence of the symmetry of the 
graph and its automorphism group. 

Symmetries of a graph, given by its automorphism group, can be inherited by the 
evolution operator. Using the irreducible representations of the automorphism group, we 
derive conditions such that quantum walks defined on this graph must have infinite hit- 
ting times for some initial states. Symmetry can cause the walk to also be confined to a 
subspace of the original Hilbert space for cartain initial states. We show that a quantum 
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walk confined to the subspace corresponding to this sjmimetry group can be seen as a 
different quantum walk on a smaller quotient graph. We give an explicit construction of 
the quotient graph for any subgroup H of the automorphism group. The automorphisms 
of the quotient graph which are inherited from the original graph are the original auto- 
morphism group modulo the subgroup H used to construct it. We conjecture that the 
existence of a small quotient graph with finite hitting times is necessary for a walk to ex- 
hibit a quantum speed-up. Finally, we use symmetry and the theory of decoherence-free 
subspaces to determine when the subspace of the quotient graph is a decoherence-free 
subspace of the dynamics. 
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Chapter 1 



Introduction and preliminaries 
1.1 Quantum computing 

In the course of the 20th century, quantum mechanics was estabhshed as a fundamental 
theory of physics. Every physical theory is built (or needs to be built) within the frame- 
work of quantum mechanics. But quantum mechanics can seem very counter-intuitive 
even to an expert. Superposition, interference and entanglement are some of the main 
and bizarre aspects of quantum mechanics. These aspects can lead to many spectacular 
and often puzzling effects in physical systems. One of the main goals of the field of quan- 
tum computing and quantum information is to use these non-classical aspects to build 
more powerful computers, communication devices and cryptographic systems. 

The invention of the transistor by Bardeen, Brattain and Shockley has revolutionized 
computer hardware. Computers have since been becoming more powerful as more and 
more transistors are made to fit on an integrated circuit (IC). In 1965, Gordon Moore 
characterized this growth in his famous Moore's law which states that the number of 
transistors on an IC doubles in number roughly every two years. This means that the 
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power of the computer doubles for the same cost roughly every two years. This law has 
held sway for many years and seems to be valid even today. However, experts predict 
that its rule will have run its time sometime in the next twenty years. As electronic 
devices become smaller they will eventually run into fundamental difficulties of size due 
to quantum effects. Conventional fabrication techniques are going to have a problem 
when these quantum effects come into play. New technologies will have to replace them 
to make progress. Quantum computing is a paradigm which can provide a solution to 
this problem. This is another reason why research in quantum computing and quantum 
information is critical. 

The original idea of a quantum computer can be traced back to simulation of quantum 
systems. Any physical system adheres to the laws of quantum mechanics and hence it 
is inherently a quantum system. Thus, understanding quantum systems is essential to 
physics. But any reasonable and large-scale model of such a system is hard to simulate on 
a classical computer since it requires exponential resources. In the 1980s, it was realized 
that an appropriate system which might efficiently simulate any quantum system is one 
which is based on quantum mechanics. The essential idea of a such a quantum computer 
was suggested by Manin [56], Benioff [11] and Feynman [31]. This idea was taken to the 
next level by Deutsch, who introduced the notion of a quantum computer as a universal 
computing machine [23, 24]. An early application of these concepts was in quantum 
cryptography [12], where a cryptographic protocol was proposed which has since been 
shown to be unconditionally secure [62]. The first instance of a computational problem 
that could be solved faster by a quantum computer than by a classical computer was 
given by Deutsch and Josza [25] . The next problem which has a faster quantum algorithm 
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than a classical one, was given by Simon [75]. Though these algorithms demonstrated 
the power of quantum computing, the problems they solve are rather artificial in nature. 
The discovery of quantum algorithms for prime factoring and discrete-logarithm by Shor 
[74] was the first example where quantum algorithms work faster than the best known 
classical algorithms for problems of great practical importance. Since then there has 
been a surge of interest in quantum computing and quantum information which led to 
the development of more areas, such as quantum Shannon theory, which was built from 
the original ideas of quantum teleportation and superdense coding [13] and quantum error 
correction. 

One of the most important goals of quantum computing is the design of fast algorithms 
for computational problems. A quantum algorithm for database search which works 
faster than any classical algorithm was given by Grover [35]. The algorithms of Grover 
[35] and Shor [74] are among the famous examples of quantum algorithms which have a 
speed-up over classical algorithms. These two algorithms are very different in structure: 
Grover's algorithm exploits an invariant two dimensional subspace within the search 
space, while Shor's algorithm exploits the properties of the quantum Fourier transform 
(QFT). Efficient algorithms for a class of problems called the hidden subgroup problem 
(HSP) (factoring and discrete-logarithm belong to this class) use the QFT [53]. The QFT 
is useful for the Abelian version of the HSP and for some non- Abelian groups [68] . But for 
other non- Abelian groups, most notably, the symmetric group and the associated HSP- 
the graph isomorhism problem, the power of the QFT seems to be limited. Grover's 
algorithm, although very useful in many search problems, gives only a quadratic speed 
up, and hence a straightforward application of this algorithm is not very efficient for the 
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HSP. This is because it ignores structure in the problem which can be used to obtain a 
speed up. New concepts and tools might be needed to design algorithms to solve these 
problems. In this context, quantum walks might provide such tools. A goal of this thesis 
is to explore the properties of quantum walks which might be useful in designing fast 
algorithms. 

We begin by giving a brief introduction to the basic concepts and definitions of quan- 
tum computing which are used in this thesis. To get a more detailed explanation refer to 
the book by Nielsen and Chuang [61]. 



1.2 The rules of quantum mechanics 

State space. An isolated quantum system can be described by a quantum state which is 
a vector in a complex vector space endowed with an inner product (i.e., a Hilbert space). 
Given an orthonormal basis {|^)} i G S, for this Hilbert space, the state of the system 
\a) can be written as \a) = 5^jCii|i), such that X^JciiP = 1. {a\ denotes the complex 
conjugate of the state and (a|/9) denotes the inner product between two states \a) and 

Unitary evolution. The evolution of a closed system can be described by a unitary 
transformation. In other words, the state of the system at two different times can be 
related by a unitary operator. If \a) and \j3) are the states at times ti and t2 respectively, 
then 

\P) = U\a), (1.1) 

4 



where U is unitary i.e., UU^ = I, W is the transpose complex conjugate of U. Note that 
the norm of the state does not change after the evolution, i.e., |(/3|/3)| = |(a|i7W|Q;)| = 



a a) ■ 



Measurement. Generalized quantum measurements are described by a set of mea- 
surement operators {Mj} which can act on the Hilbert space of the system and satisfy the 
completeness relation J2i M^Mi = I. The index i denotes the measurement outcome and 
when a measurement is performed any one of these outcomes can occur with a certain 
probability. Given a state of the system, say |q;) before the measurement, an outcome k 
occurs with a probability {a\MlMj^\a) and if this outcome k occurs, then the resulting 
state of the system is, 

fel / ^^ 2) 



'{a\MlMk\a) 

Projective measurements are a special case of generalized measurements in which the 
measurement operators are projectors. An operator P is a projector if = P and 
P = P^. Thus, in such measurements, the completeness relation becomes X^jPj = I, 
where Pj are the measurement operators. The operators Pj also satisfy orthogonality 
relations PjPj = Sij. 

Density operators. Suppose the exact quantum state is not known, but rather it 
is known that the state belongs to a set where i indexes the set. If each could 
be the state with a probability pi, then the density operator (also known as the density 
matrix) for the system is defined as 



i 



(1.3) 
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If the state of a quantum system is known exactly (1^/;), say), it is called a pure state. The 
density operator corresponding to a pure state is p = If the state is not pure, 

it is called mixed. The unitary evolution and measurement postulates can be written in 
the density operator picture in the following way. Given that the density matrices of a 
closed quantum system at times ti and t2 are pi and p2 respectively, this evolution can 
be related by a unitary operator as 

P2 = UpiU^. (1.4) 

Now, given a measurement described by measurement operators {Mi} (such that J2i ^l^i = 
I) , the outcome k occurs with a probability p{k) = Ti{MlMkp) and the state of the system 
is then Mi;pMI/Ti{mIMi;p). For a projective measurement with operators {Pi}, an out- 
come k occurs with probability p{k) = Tr(Pfep) and the resulting state is PkpPk/T^{PkP)- 
Density operators are characterized by the following two properties. 

1. The trace of a density matrix is unity. 

2. A density matrix is a positive matrix. 

1.3 Outline 

The rest of this thesis is organized as follows. In Chapter 2, we define classical random 
walks and then give the definitions of quantum walks. Then we define Cayley graphs and 
give a detailed description of the hypercube as an example of a Cayley graph. 
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In Chapter 3, we define a notion of hitting time and provide a formula for the hitting 
time of a quantum walk on a graph. We then give the results of simulations for various 
cases with and without decoherence for the hitting time on the hypercube. We also 
provide evidence of infinite hitting times and then derive rigorously the conditions under 
which a quantum walk has infinite hitting times. 

In Chapter 4, we relate the notion of infinite hitting times to symmetry of the graph 
and give examples of graphs with sufficient symmetry to have infinite hitting times. We 
also provide the necessary definitions and results in the representation theory of finite 
groups used in the analysis. 

In Chapter 5, we use symmetry again to explain fast hitting times. We develop 
the notion of a quotient graph and use it to show that quantum walks on a graph with 
symmetry are actually walks on a smaller quotient graph. Quantum walks with symmetry 
have invariant subspaccs which can be used for algorithm design. 

In Chapter 6, we make use of symmetry to find conditions on the decoherence such 
that the subspace of a quotient graph will lie in a decoherence-free subspaces of the 
dynamics. 
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Chapter 2 



Random and quantum walks on graphs 



Quantum walks were formulated in studies involving the dynamics of quantum diffusion 
[30], but the analysis of quantum walks for use in quantum algorithms was first done by 
Farhi and Gutmann [27] . The motivation for the development of quantum walks is three- 
fold. First, algorithm design in classical computer science has benefitted enormously with 
the advent of randomized algorithms. The best known algorithms for some important 
problems like 3-S AT (3-satisfiability) are based on random walks [59] . A natural question 
is whether quantum walks would perform better for these problems. Second, as mentioned 
above existing quantum algorithms for the so called hidden subgroup problem, based on 
the quantum Fourier transform, like Shor's algorithm and related algorithms [74], do not 
seem to be effective for certain non-Abelian problems like the graph isomorphism problem. 
Therefore, there is a need for a new class of algorithms to tackle these problems, and 
quantum walk based algorithms may provide a new approach. Finally, a large number of 
problems in computer science can be reformulated as graph related problems and having 
a quantum formalism which is specific for graphs may be useful in designing algorithms 
for these problems. 
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Quantum walks have been already been used in quantum algorithms which have a 
speed-up over the corresponding classical algorithms. Childs et al. showed in [19] that 
a continuous quantum walk would find a certain final node exponentially faster than a 
classical random walk on the so-called "glued trees" graph. The classical walk tends to 
take exponentially longer due to the large number of vertices in the middle of the graph, 
where as a quantum walk would traverse it in a superposition of paths and move toward 
the final vertex much faster. Shenvi et al. [72] showed that a quantum walk on an 
unsorted database (represented as a hypercube) has a quadratic speed up over a classical 
algorithm. The algorithm begins in a superposition of all the vertices and coin states and 
proceeds by applying the "Grover" coin at every vertex except the final vertex where it 
applies another coin (— /, which they point out in the paper, is an arbitrary choice since 
other coins worked just as well in numerical simulations.) This makes the walk converge 
to the final vertex in 0(\/iV) steps, where N is the number of vertices. Quantum search 
for a marked item in N items arranged on a d-dimensional N^/''' x N^/''' x . . . N^/'^ grid 
was analyzed using a continuous-time walk by Childs and Goldstone in [20] and using 
a discrete-time walk by Ambainis et al. in [7]. This is one of the few places where 
discrete walks are better. An algorithm based on a continuous-time walk worked in 
0{VN) time steps for d > 5, in 0{y/N log N) for d = 4 and had no speed up in d = 2, 3 
whereas a discrete walk took 0{-\/N) for d > 3 and 0{-\/N\ogN) for d = 2. Ambainis 
[4] has applied quantum walks to the element distinctness problem i.e., the problem of 
determining whether A'^ elements in a given set are all distinct or not. This algorithm is 
based on a discrete-time quantum walk on the Johnson graph. Other algorithms based 
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on quantum walks are an algorithm for matrix product verification [17], triangle finding 
[54] and group commutativity testing [55]. 

Several quantities of interest have been defined for quantum walks analogous to clas- 
sical walks in [2], such as mixing time, sampling time, filling time and dispersion time. In 
[2], a lower bound on the mixing time of a discrete quantum walk on the A^-cycle is found 
and shown to be at most polynomially faster than that of the classical walk. Hitting time 
is another important quantity for classical walks on graphs. Two definitions of hitting 
time are given in [38] and an upper bound for one of them was found for the walk on 
a hypercube. A different definition of hitting time is given in [47], where the unitary 
evolution of the discrete walk is replaced by a measured walk. In such a walk, after the 
application of the unitary evolution operator, a measurement is performed to see if the 
particle is in the final vertex or not. In [48] the phenomenon of infinite hitting times is 
analyzed and it was shown that graphs with sufficient symmetry can have infinite hitting 
times for certain initial states. This is a purely quantum phenomenon and does not have 
a classical analogue. Cayley graphs on the symmetric group are shown to be examples 
of graphs which have this symmetry and hence have an infinite hitting time for certain 
starting states. The theory of irreducible representations is used to estimate the amount 
of degeneracy that a given group of symmetries produces. The use of representation 
theory to explain aspects of quantum walks on certain classes of graphs was also done 
in [33], where the behavior of mixing times of Cayley graphs on the symmetric group is 
explained based on its irreducible representations. 

Quantum walks have applications other than in the design of new algorithms. In 
[21, 22], quantum walks on weighted graphs have been used to efficiently transfer quantum 
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states with perfect fidelity. Here and in [28], symmetry has been used to demonstrate a 
class of graphs on which a continuous time quantum walk reduces to the walk on a line 
(quantum wire) like in the case of the glued-trees graph. 

Quantum walks come in two distinct flavors: discrete-time and continuous-time. The 
main difference between them is that discrete time walks require a "coin" — which is just 
any unitary matrix — plus an extra Hilbert space on which the coin acts, while continuous 
time walks do not need this extra Hilbert space. Aside from this, these two versions 
are similar to their classical counterparts. Discrete-time quantum walks evolve by the 
application of a unitary evolution operator at discrete time intervals, and continuous- 
time walks evolve under a (usually time-independent) Hamiltonian. Unlike the classical 
case, the extra Hilbert space for discrete-time quantum walks means that one cannot 
obtain the continuous quantum walk from the discrete walk by taking a limit as the time 
step goes to zero. Although there is no natural limit to go from the discrete to continuous 
walks for general graphs, for the quantum walk on the line [77] offers a treatment of this 
limit, where it is possible to meaningfully extract the continuous-time walk as a limit of 
the discrete-time walk. The dynamics of quantum walks of both types has been studied in 
detail for walks on an infinite line — for the continuous-time case in Refs. [19, 27, 18, 45, 46] 
and for the discrete-time case in [60, 10, 14, 15, 16]. There has also been considerable 
work on other regular graphs. The A'^-cycle is treated in [2, 80], and the hypercube in 
[72, 57, 38, 47, 48, 49, 50]. Quantum walks on general undirected graphs are defined in 
[42, 5], and on directed graphs in [58]. Reviews of quantum walks include an introductory 
review by Kempe in [37], and a review from the perspective of algorithms by Ambainis 
in [5]. The role of symmetry in quantum walks has been analyzed in [47, 48, 49]. It has 
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been shown in [48] that the evolution operator of a quantum walk inherits symmetries 
from the automorphisms of the graph and this leads to degeneracy in the operator. This 
degeneracy can determine subspaces to which the walk remains confined. Thus, the walk 
never explores some parts of the Hilbert space, leading to the phenomenon of infinite 
hitting times. In [49], it was shown that due to the symmetry of the graph, certain 
subspaces which confine the walk have fast hitting times, and moreover these subspaces 
are the corresponding Hilbert spaces of a different graph - a quotient graph. Thus, the 
walk is confined to smaller graph due to symmetry. The exponential speed-up observed 
for quantum walks in [19, 72] on the "glued-trees" graph and the hypercube respectively 
is due to the fact that the quantum walk is on a smaller quotient graph. 

Decoherence in quantum systems can be broadly defined as any process that destroys 
quantum coherence. Aharanov [1] first considered quantum walks with measurements. 
Decoherence was later considered in quantum walks mainly to demonstrate the classical 
behavior of such decohering quantum walks. Quantum to classical transition due to 
decoherence has been considered for the walk on the line in Refs. [14, 15] and for the 
hypercube in [3, 41]. Kendon et al [43] provide a treatment of this transition for any 
quantum walk. More recently, decoherence has been shown to be beneficial for fast 
mixing behavior of quantum walks. For the continuous-time walk on the hypercube, 
decohering quantum walks have been shown to mix faster [63, 64]. Kendon [44] gives a 
review of the work done in this field so far, focusing mainly on decoherence. Studies of 
decoherence in quantum walks is useful not only to understand the quantum behavior of 
the walk but also for possible implementations. Any implementation of a quantum walk 
will involve having to deal with decoherence. Such studies may be useful in determining 
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the kind of decoherence that has httle or no effect on the useful properties of the walk. 
Some implementation schemes proposed so far are in Refs. [26, 29, 43, 70, 69]. 



2.1 Random walks 

A simple random walk on an undirected graph is defined as the repeated application of a 
stochastic matrix P, where P{i,j) = l/diifi and j are connected and otherwise, where 
P{i, j) is the probability to go from vertex i to vertex j and di and dj are the degrees 
of these vertices. Such a process is called a Markov chain. If the graph is connected 
and non-bipartite then a fundamental property of Markov chains is that the distribution 
tends to a stationary distribution (tt) which is independent of the initial distribution 
If the graph is regular (if every vertex is connected to the same number of other vertices), 
then this final distribution is uniform over all the vertices. By contrast, quantum walks 
do not converge to a final distribution since norms of states do not change under a 
unitary operation and hence the distance between the states describing the system does 
not converge to zero ([2], [37]). The induced probability distribution does not converge 
either, but it turns out that the time-averaged probability distribution converges. The 
time averaged distribution can be defined as 



where |ao) is the initial state, Pt is the instantaneous distribution and v is any vertex of 
the graph. 




T-l 



(2.1) 



t=o 
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The rate of convergence to this final distribution can expressed in terms of many 
quantities, but the one used commonly is mixing time. Mixing time is defined as 

Me = min{r|Vi > T,/ : - 7r|| < e}, (2.2) 

where | |p — q| | = \Pi — Qi\ is the total variation distance between the distributions p and 
q. Mixing time for a classical walk is related to the spectral gap (the difference between 
the largest and the second largest eigenvalues of P) in the following way: 

(A.-A;)log2. S ^ (a7^('"."'°'5''" +log^-'). (2.3) 

where Ai is the largest eigenvalue and A2 is the second largest eigenvalue. If the graph 
is regular, it turns out that the largest eigenvalue of P (i.e., Ai) is 1. The above relation 
connects mixing time and the second largest eigenvalue of the stochastic matrix. The 
classical mixing time has been found for a number of graphs. For the A?"-cycle the mixing 
time of a simple random walk can be calculated to be 0(iV^ log(l/e)). The mixing time 
for the hypercube turns out to be 0(dlog (ilog(l/e)), where d is the dimension of the 
hypercube. The hypercube is an example of a class of graphs called expander graphs for 
which it has been shown that a classical random walk is rapidly mixing. But quantum 
walks on the other hand do not have good mixing properties on the hypercube. The 
discrete-time quantum walk has a mixing time of at least 0((i^/^/e) [57], where d again 
is the dimension of the hypercube. But a quantity that looks promising for quantum 
walks on the hypercube is hitting time. Hitting time measures the average time it takes 
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for a walk to reach a certain vertex from a given starting vertex. The analysis of hitting 
times under different situations is one of the aspects of this thesis. In order to motivate 
the analysis of hitting times of quantum walks on the hypercube, we briefly present the 
following classical random walk based algorithm. 

Classical random walks have been used in many randomized algorithms in computer 
science. One of the applications of random walks to algorithms is for the 3-satisfiability (3- 
SAT) problem where the best known classical algorithm is based on the hitting time of a 
random walk on the hypercube. The 3-SAT problem can be defined as follows. Consider a 
set of n literals {xi,X2, ■ ■ ■ , Xn}, each taking a value of or 1, and a set of clauses such that 
each clause contains a logical OR of only three literals (or their negations). The problem 
of 3-satisfiability consists of finding an assignment for the literals such that each of the 
clauses is satisfied (i.e., each has a value 1). For example, (a;iV-iX2V-ia;3) A(-ixiVa;2Va;4) 
is an expression with two clauses, where A and V denote the binary operations of AND 
and OR and -^x denotes the negation of x. (xi, X2, X3, 2:4) = (1,0,0, 1) is an assignment 
of the literals that satisfies both the clauses and thus represents a solution. The classical 
random-walk based algorithm takes 0(1.329..."poly(n)) steps [65, 73]. It consists of the 
following steps. 

1. Choose a random initial assignment for the literals. 

2. Repeat 3n times 

(a) If all the clauses are satisfied, then stop. 

(b) If not, pick one unsatisfied clause, choose one literal uniformly at random and 
flip it. 
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This algorithm is a random walk on the hypercube since an initial assignment of the 
literals represents a vertex on the hypercube and flipping a bit of this vertex represents 
moving along an edge connected to it. Thus analyzing the behavior of hitting times 
of quantum walks on the hypercube may give us important clues to design quantum 
algorithms for the 3-SAT. However, in the above algorithm, the walk is on a directed 
hypercube since after moving along an edge to a new vertex, the literal that was flipped 
need not occur in an unsatisfied clause. This means that it cannot be flipped at the new 
vertex and therefore, there is no path from this vertex to the original vertex. This directed 
nature of the graph causes problems when we try to make this a quantum algorithm. A 
quantum walk as with any quantum computation, needs to be unitary. It has been 
shown in [58] that on a directed graph, it is possible to define a unitary walk if and only 
if the graph is reversible. A directed edge of a graph going from vertex Vi to Vj is called 
reversible if there is a path from Vj to Vi. A graph is called reversible if every edge in it is 
reversible. Thus, given an expression for 3-SAT, one could define a unitary quantum walk 
if the corresponding directed graph of the problem is reversible. This may be overcome 
by defining a non-unitary walk, one that involves measurements, if it is not reversible. 
But the next step, which is to design an algorithm based on this walk, is not clear. There 
does not exist a quantum walk algorithm for 3-SAT presently, but a quantum algorithm 
based on Grover search and amplitude amplification which works in 0(1.153..."poly(n)) 
steps is presented in [8]. 
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2.2 Quantum walks— discrete and continuous 



Quantum walks as noted earlier, are primarily of two types. Depending on the way 
the evolution operator is defined, they can be either discrete-time or continuous-time 
quantum walks. These two definitions of quantum walks are not exactly equivalent to 
the two types of classical random walks. While they are both based on the classical 
definitions, unlike the classical case the discrete quantum walk does not reduce to the 
continuous walk when we let the time step between repeated applications of the unitary 
tend to zero. This is because discrete-time walks need an extra Hilbert space, called the 
"coin" space (from the idea that one flips a coin at each step to determine which way 
to walk), and taking the limit where the time step goes to zero does not eliminate this 
Hilbert space. Therefore, the properties of discrete and continuous walks are different. 
Though there is no obvious reason why one should be preferred, in some cases it has been 
shown that coins make these walks faster [7] . 

2.2.1 Discrete-time walks 

A discrete-time quantum walk can broadly be defined as the repeated application of 
a unitary evolution operator on a Hilbert space whose size depends on the graph. This 
Hilbert space usually consists of the space of possible positions (i.e., the vertices) together 
with the space of possible directions in which the particle can move along from each vertex 
(the coin space). All the concepts which we develop in this thesis such as a formula for 
hitting time of a dicrete-time quantum walk, infinite hitting times, effect of symmetry 
and quotient graphs are applicable to any undirected graph. However, in most of our 
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simulations we consider d-rcgular, undirected and cZ-colorable graphs. This is because the 
structure of quantum walks reduces to something more manageable in this case. We will 
briefly review the definitions of these graph properties. 

A regular graph is one where every vertex is connected to the same number d of other 
vertices. This number is called the the degree of the graph. A graph is undirected if for 
every edge between vertices A and B going from A to B, an edge goes from B to A as well. 
In this case, we identify the edge from A to B with the edge from B to A, and consider 
them a single edge. A regular, undirected graph with vertices of degree d is considered 
d- colorable if the edges incident on every vertex can be numbered 1 through d such that 
every edge between two vertices has the same number at either end. Not all d-regular 
and undirected graphs can be d-colored. A simple example is the triangle graph where 
N = 3 and d = 2. (See Fig. 2.1.) For o?-regular, undirected and c?-colored graphs, the 
Hilbert space of the walk is W^T-C^, i.e., the tensor product of the position and direction 
(or coin) space. The evolution operator U is given hy U = S{I ® C), where S is called 
the shift matrix and C is the coin matrix. The shift matrix encodes the structure of the 
graph and is very similar to its adjacency matrix. The vertices, numbered |0) through 
|iV — 1), are basis states for the vertex Hilbert space TiP and the set of all directions from 
each vertex, numbered |1) through \d), are basis states for the coin Hilbert space H"^. In 
this basis, the shift matrix for the graph can be given the explicit form: 

V i 

where v{i) is the vertex connected to v along the edge numbered i. 
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The coin matrix C acts only on the coin space, and "flips" the directions before the 
shift matrix is applied. Then S moves the particle from its present vertex to the vertex 
connected to it along the edge indicated by the coin direction. Though C can be any 
unitary matrix, usually coins with some structure are considered. The coins that we used 
in our previous analysis are the Grover coin G and the Discrete Fourier Transform (DFT) 
coin D. The matrices for these coins are given by: 



G = 2\^){^\ - I 



2 _ 1 2 
d ^ d 



V i 



- 1 



(2.4) 



and 



D 



1 1 

1 u; 



to 



\1 u'^-^ 

where |^) = ^ |i) and u) = exp{2m / d) . 



1 

d-l 



(2.5) 



2.2.2 Decoherence in discrete-time walks 

Decoherence in any quantum system is a process which destroys quantum superpositions. 
Any quantum system will have decoherence (quantum noise). Thus any practical imple- 
mentation scheme of quantum walks must deal with decoherence. In general, it can be 
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Figure 2.1: Examples of 2-colorable (square) and non-2-colorable (triangle) regular graphs 
of degree 2, where the colors are numbered 1 and 2. 

thought of as any completely positive and trace preserving (CPTP) map. An arbitrary 
CPTP map can be represented using the operator sum (OSR) or Kraus representation. 

P(p) = ^i,pil, (2.6) 

i 

where the operators Ai are the Kraus operators and satisfy the relation: ^1^? = I- 
Thus in a discrete-time quantum walk with decoherence, we can assume that a step of 
the evolution is of the form 

£{p)=VoU{p), (2.7) 

where U{p) = UpW as before. 
2.2.3 Continuous-time walks 

Continuous time quantum walks were defined by Farhi and Gutmann in [27]. For an 
undirected graph G{V,E), the unitary evolution operator is defined as U = exp(iHt), 
where H is obtained from the adjacency matrix of the graph. Here again, the vertices of 

20 



the graph form a basis for the Hilbert space on which U is defined. This gives rise to the 



foUowing Schrodinger equation: 



(2.8) 



This walk has a structure very similar to that of continuous time Markov chains. H is 
defined as 



where 7 is the jumping rate from a vertex to its neighbor i.e., the transitions between 
connected vertices happen with a probability 7 per unit time. 

But for a regular graph we can take H to be the adjacency matrix because di = d, 
where d is the degree of the graph. This means that the Hamiltonian can be written as 
H = ^{D — A), where D = dl and A is the adjacency matrix of the graph. The matrix 
D would lead to a trivial phase factor and can be dropped. If is a symmetric matrix 
(and hence IJ unitary) if the graph is undirected. Therefore, for a regular and undirected 
graph, the adjacency matrix H, which acts as the Hamiltonian, is of the form: 



~7 ^ 7^ J if nodes i and j connected 



i / j if nodes i and j not connected 



(2.9) 



da i = j 




1 Hi and j share an edge. 



(2.10) 



otherwise. 
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As can be seen, this walk has no coin and so the Hilbert space on which U acts is only 
the vertex space HP. 

2.2.4 Decoherence in continuous-time walks 

Continuous-time quantum walks evolve by an application of a Hamiltonian as defined in 
Eq. (2.9). Decoherence in this scenario would give rise to an evolution which can be 
described using the Lindblad semigroup master equation as 

p{t) = [H, Pit)] + Y,iL^pit)Ll + LlUpit) + pmlu), (2.11) 

i 

where Lj are the Lindblad operators. The Lindblad operators represent the environmental 
interactions which lead to decoherence. 

2.3 Cay ley graphs 





011 


/ 




U1U 




100 






c 


/ 

00 


[/ 

001 


101 



Figure 2.2: The hypercube in three dimensions as an example of a Cayley graph. 

Cayley graphs are defined in terms of a group G and a set S consisting of elements from 
G such that the identity element e ^ S. Given G and S, the resulting (right)-Cayley graph 
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r(G, S) is one whose vertices are labeled by the group elements, i.e., there is one vertex 
for every group element, and two vertices g and h are connected by a directed edge from g 
to h if g~^h G S, (see [34]). Another way to look at this definition is that from any vertex 
5 of a Cayley graph, there are l^l outgoing edges, one to each of the vertices gs, Vs G 5. A 
Cayley graph will be connected if and only if the set S is a generating set for G, and it will 
be undirected if G -S, Vs G S. The degree of such a graph is \S\, the cardinality of the 
generating set. Finally, a d-regular Cayley graph can be d-colored if = 1, Vs G 5, i.e., 
= s. Examples of Cayley graphs on which quantum walks have been studied include 
the line T{Z, {1, —1}); the cycle r(Z„, {1, — 1}); the hypercube T{Z2,X) where the set 
X is the set of canonical generators {(1, 0, 0, • • • , 0), (0, 1, 0, • • • , 0), . . . , (0, 0, 0, • • • , 1)}; 
and the graph on the symmetric group r(S'„, Y), where y is a generating set for Sn- Let 
us look at the hypercube as an example of a Cayley graph where quantum walks have 
been extensively studied (see Fig. (2.2)). The hypercube has 1-2^2 I = 2" vertices each with 
a degree of \X\ = n. The vertices can be labeled by an n-bit string from (0,0,- • • ,0) 
through (1,1,-- - ,1). Two vertices are adjacent if they differ only by a single bit. Vertex 
V is connected to n vertices given hy v®s, ^se X, where © stands for the bit- wise XOR 
of the bit strings v and s. One important property of the hypercube is that it can be 
n-colored, since s s = e, Vs G X where e = (0, 0, • • • ,0) is the identity element. The 
unitary evolution operator for a discrete walk on the hypercube becomes J7 = S{I C), 
where S has the form 

S = ^^\v® s){v\ (g) \s){s\. 

S V 
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Since the vertices of the hypercube are bit strings, and adjacent vertices are those that 
differ by one bit, the shift matrix of the discrete walk on the hypercube has a natural 
form given by 

S = X(g)/(g)---(g)/(g) +/(g)X(g)---(g)i(g)|s2)(s2| 

+ ... + 7(g)7(g)---(g)X(8)|s„)(s„|, (2.12) 

where X stands for the Pauh ax operator. This structure of S reflects the property of the 
hypercube that moving along an edge from v corresponds to flipping one bit of v. This 
structure is also useful in determining its group of symmetries as we shall see later. 
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Chapter 3 



Hitting times 

3.1 Classical hitting time 

Given a regular undirected graph and a particle which starts at some vertex, the classical 
random walk is defined as before. At each vertex, the particle moves along any edge 
incident on the vertex with some predefined probability. This procedure is then repeated 
at the new vertex. The walk continues until the particle arrives at ("hits") a certain 
vertex (called the "final vertex") for the first time. The hitting time is defined as the 
average time until the particle hits the final vertex: 



where t{v) is the hitting time given that the walk starts at vertex v and Pv{t) is the 
probability that the particle hits the final vertex for the first time at time step t (first 
crossing probability) given that it was at u at t = 0. 

Let us now specialize to the case of the hypercube, where the the final vertex is 
assumed to be 11 • • • 1. We would like to find the hitting time starting from 00 - • • 0. 



oo 




(3.1) 



t=0 
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For the classical walk on the hypercube, one can arrive at a recursive relation involving 
the hitting time. First, from the symmetry of the hypercube one can conclude that the 
hitting time depends only on the Hamming weight of the starting vertex rather than the 
vertex itself. The Hamming weight is the number of I's in the string of bits. At Hamming 
weight X, there are = nl/xl{n — x)\ vertices. The probability to walk to a vertex with 
weight X + 1 is {n — x)/n, and the probability to walk to a vertex with weight x — 1 is 
x/n. So, if t{x) denotes the hitting time starting at any vertex with Hamming weight x, 
then 

x) = r x + 1 +1, 3.2 

n n 

with the boundary condition r(ra) = 0. This simplifies to 

where A(x) = r(x) — t{x + 1). Using this recursive formula, we obtain 

n— 1 n—1 v^x— 1 i i 

x=0 x=0 

This sum can readily be evaluated for reasonable sizes of n and in fact this sum scales 
as 2". We use this expression to compare the classical hitting time to the quantum 
hitting time. We define the hitting time of a quantum walk next. 
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3.2 Hitting time for quantum walks 



The hitting time Th of a classical random walk is defined as the average time for the 
walk to hit a designated 'final' vertex Vf given that the walk began with some initial 
distribution pf. 



where p{t) is the probability of being in the final vertex for the first time at time step t. 
In order to carry this notion of hitting time over to the quantum case, we need to make 
the meaning of pit) more precise. In particular, we need to define clearly what "for the 
first time" means for a quantum walk. As described in [47], we do this by performing 
a projective measurement of the particle at every step of the walk to see if the particle 
has reached the final vertex or not. The measurement M which is used has projectors P/ 
and Qf = I — Pf representing the particle being found or not found at the final vertex, 
respectively. The projector is defined Pf = \xf){xf\ (g) Ic, where \xf) is the final vertex 
state and Ic is the identity operator on the coin space. Using this definition, each step of 
the measured walk consists of an application of the unitary evolution operator U followed 
by the measurement M. 

Now we can use the same expression (3.5) for the hitting time, where the probability 
p{t) becomes 



oo 




(3.5) 



t=o 



Pit) = T^^iPfUiQfUf-'poiU^Qft'U^Pf}. 



(3.6) 
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To explicitly sum the series in Eq. (3.5) using tiie expression for p{t) in Eq. (3.6), we 
rewrite the expression in terms of superoperators (linear transformations on operators) M 
and y, defined by 



Afp = QfUpWQf 

yp = PfUpWPf. (3.7) 

In terms of M and y, p{t) = TT{yM^~^ po} ■ We introduce a new superoperator 0{l) 
which depends on a real parameter I: 

oo 

o{i) = ij2m'-\ (3.8) 

t=i 

which is a function of the parameter I. The hitting time now becomes 



Th = j^Tr{yO{l)po} 



(3.9) 



If the superoperator I — IM is invertible, then we can replace the sum Eq. (3.8) with 
the closed form 

0{l) = 1{I - lM)-\ (3.10) 

(The case when I — lj\f is not invertible is discussed in detail later.) The derivative in 
Eq. (3.9) is 

^(1) = {I-J\f)-'+J\f{l-J\f)-^ = (I-M)-^. (3.11) 
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This gives us the following expression for the hitting time: 



(3.12) 



To evaluate Eq. (3.12), we write these superoperators as matrices using Roth's lemma 
[66]. As shown in [47], we can then vectorize the density operators and operators on states, 
and write the action of superoperators as simple matrix multiplication. Any matrix can 
be vectorized by turning its rows into columns and stacking them up one by one, so that 
Si, D X D matrix becomes a column vector of size D^. For example: 

an 
ai2 
ai3 

021 
022 
023 
031 
032 
y033 J 

Consequently the superoperators become matrices of size D"^ x D"^. This method of 
vectorization takes operators on one Hilbert space H to vectors in another Hilbert space 
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012 


Ol3 


021 


022 


023 


^031 


032 


O33 j 



H' = H'S>H* and so superoperators in H are operators on H'. Note that a basis {|itij)} 
for H' can be obtained from a basis {\vi)} for H by defining 




Vi) (g) \Vj)*. 



(3.13) 



For our superoperators J\f and y we then get 



{^fpr = 



{yp) 



V 



{PfU)®{Pfur p\ 



(3.14) 



Let N = {QflJ) ® (QfU)* and Y = {PfU) ® (P/f/)*. The hitting time becomes 



Using this vectorization transformation, we treat the superoperators as operators on a 
larger Hilbert space and thus can find their inverses (when they exist). But, the expression 
in Eq. (3.15) is not always well defined because the matrix I — N may not be invertible. 
We will show that when it is not invertible, it means that the hitting time becomes infinite 
for some initial states, and vice versa. This property of quantum walks having infinite 
hitting times does not have a classical analogue. 



T;, = r-(Y(I-N)-V)- 



(3.15) 
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3.3 Other definitions of hitting time 

Two definitions of the quantum hitting time were given in [38] : one-shot hitting time and 
concurrent hitting time. The one-shot hitting time is defined for an unmeasured walk. It 
is the time at which the probability of being in the final state is greater than some given 
value. More precisely, given some probability p, the one-shot hitting time is defined as 
the lowest time rsh(p) such that 

\{xf\U^-^^P^\xo)\'^ >p (3.16) 

where Xf and xq are the final and initial states, and U is the evolution operator (as defined 
above). Essentially this same definition of hitting time was used in the analysis of the 
continuous-time walk on the hypercube in [3] . This definition is useful if it is known that 
at some time the probability to be in the final state will be higher than some reasonable 
value; but for a general graph, this is not guaranteed. 

The concurrent hitting time, by contrast, is defined for a measured quantum walk. 
Given a probability p, it is the time Tc{p) such that the measured walk has a probability 
greater than p of stopping at a time less than rdp). It has been proved that the con- 
current hitting time is Tc{p) = | for p = ^{j^j^^) for a hypercube of dimension n (for 
the symmetric initial condition used in this paper) . Since we consider only the measured 
quantum walk in our analysis, we compare our numerical results to the numerical sim- 
ulation of the concurrent hitting time and the bound on it derived in [38]. In the next 
section, this is redefined in terms of the residual probability 1—p and plotted against the 
hitting time defined in the previous section. 
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If we think of quantum walks as a possible route to new algorithms, then the con- 
current hitting time corresponds to the time needed to find a solution with probability 
greater than p. The definition of hitting time used here corresponds more to a typical 
running time for the algorithm. Both definitions could prove useful for particular pur- 
poses. In particular, both definitions give different characterizations of the distribution 
p{t) — the probability of first being in the final vertex at time t — in different ways. A 
complete understanding of the time needed to find a solution would require knowledge 
of the entire distribution p{t) (or Tc{p) for all p), which is unlikely to be achievable, in 
practice; but Th and Tc{p) (for fixed p) give different windows on this function. 



3.4 Results for the Grover coin 

We calculated the average hitting time by evaluating the above expression (3.15) in 
Matlab. Because of the multiple tensor products, the size of the matrices J\f and y 
is (2"ri)^ X (2"ri)^. As n increases, the size of these matrices increases exponentially, and 
the matrix inversion that Eq. (3.15) demands is not easy to compute. We can always 
compute an estimate of the quantum hitting time (strictly, a lower bound) , obtained by 
iterating the quantum walk for a large number of steps. To get this lower bound, we 
define rest(e) using Eq. (3.5), and summing the series up to a finite number of terms: 



Tc(l-e) 




(3.17) 



t=i 
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where the concurrent hitting time Tc{1 — e), as defined in [37], is the shortest time T for 
which 

T 

Y,pit)>l-e. 
t 

This comparison of the exact and approximate values of the average hitting time is useful, 
because the matrices in Eq. (3.15) grow faster than those used in summing Eq. (3.5), 
which means that for some graphs it may be impractical to calculate the exact value of 
T/j, but still possible to calculate the lower bound rest(e). However, for the quantum walk 
on the hypercube with the Grover coin, when the initial state of the walk is a particular 
symmetric state, it is possible to reduce the matrix sizes considerably since the walk 
remains in a lower dimensional subspace of the original space, as explained below. For 
such an initial state we can compute the exact average hitting time (Eq. (3.15)) up to a 
large number of dimensions. Figure 3.1 shows both the classical and quantum walks on 
the hypercube for dimensions up to 100. The exact average hitting time Th is plotted as a 
dotted line, and the lower bound Test(e) for e = 0.001 is plotted as a solid line. These two 
lines almost coincide in the graph, and we can see that when e is small enough summing 
the series gives a very good estimate of t^. We conjecture that this will remain true for 
more general graphs as well. In comparing classical and quantum results, the average 
hitting time for the quantum walk is a low order polynomial, whereas the classical walk 
grows exponentially with dimension, so there is a very dramatic speed-up in the quantum 
case. 

The simplification that makes this computation tractable is in the case where the 
coin-flip unitary is the Grover coin G defined in (2.4), and the specific starting state is 
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p = |0)(0| where = Y17=i K)- I* was first observed in [57] that for the walk 
on the hypercube with this initial state and the Grover coin, the state remains always in 
a 2n-dimensional subspace, where the walk is on a line with n + 1 points. This is rather 
similar to the simplification made in the classical case, when we kept track only of the 
Hamming weight of the current vertex. With this simplification, the operators in Eq. 
(3.12) reduce to {2n)^ x (2n)^ matrices obtained in [72], which makes explicit calculations 
possible even for high dimensional hypercubes. 

We write a set of basis states for this subspace as \R, 0), \L, 1), \R, 1), . . . ,\R,n — 
1), \L,n), where the first label says whether the state is "right-going" or "left-going", and 
the second label gives the Hamming weight of the state. The initial state is |-R, 0), and 
the final state is \L, n). (Note that there are no states |L, 0) or \R, n).) Restricted to this 
subspace, the S and C matrices become 

n 

S = Y^\R,x){L,x + l\ + \L,x + l) {R, x\ (3.18) 

x=0 

and 

C\L,x) = — cos LOx\L,x) + smLOx\R,x), 

C\R,x) = sina;j;|L,a;) -I- cosa;a;|i?, x), (3.19) 

where cosux = 1 — 2x/n. We see that for the walk in this subspace, the coin flip is no 
longer independent of the position; this is quite analogous to the reduction of the classical 
walk to the Hamming weight, in which the probabilities favor walking toward x = n/2. 
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How does the average hitting time compare to the concurrent hitting time for the 
walk on the hypercube? Figure 3.2 plots the estimate of the average hitting time Test(e), 
the concurrent hitting time Tc{1 — e), and the bound on Tc(1 — e) (obtained in [38]) as 
a function of e for dimensions from 10 to 20. (Both the axes are in log scale.) Figure 
3.3 plots the same for dimensions from 50 to 60. We can see from these two figures that 
both Tc{l — e) and the bound on it become less tight for higher dimensions, and both are 
much longer than the bound on the average hitting time rest(e)- This, together with the 
comparison of Tcst(e) to the exact value Th, strongly indicates that the average walk ends 
much faster than these bounds might suggest, and also that the estimate rest(e) is quite 
insensitive to the choice of e, at least for the hypercube. 
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Figure 3.1: Hitting times of Classical and Quantum walks (semi- log scale) 
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Figure 3.2: Comparison of Quantum hitting time, concurrent hitting time and the bound 
on the concurrent hitting time for dimensions 10 to 20(log scale) 

3.4.1 Hitting time in the presence of decoherence 

This decoherence can be thought of as acting on the system at every step of the measured 
walk. Thus the evolution of the walk can be decomposed into three parts~the unitary 
evolution operator followed by the decoherence (CPTP) map and then the partial mea- 
surement on the final vertex all taking place in sequence. The combined effect can be 
seen as the application of a superoperator at every step and it can be written as 



p{t + 1) 



£{p{t)) 



MoVoU{p{t)) 




(3.20) 
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Figure 3.3: Comparison of Quantum hitting time, concurrent hitting time and the bound 
on the concurrent hitting time for dimensions 50 to 60 (log scale) 



where Mj is the measurement operator with the outcome j £ {0, 1} with Mq = Pj and 
Ml = Qf- The state of the system after n steps can be written as 

p{t)=£\p{0)). (3.21) 



The formula for hitting time in Eq. (3.15) assumes a purely unitary evolution and 
partial measurement. We now modify this formula to include the effect of decoherence on 
the walk. In the presence of decoherence, the evolution without measurement is described 

by 

p{t + l) =V oU{p{t)). (3.22) 
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Therefore, we only need to replace the superoperator U to T> oU. This means that the 
vectorized quantities Y and N are modified accordingly. They become 

Yd = (P/®P;)(^ii®i*)(?7®{7*) 

i 

= {Qf®Q})iJ2^i®^t)i^(S)U*), (3.23) 

i 

and we can write D = J2i A <^ as the decoherence superoperator. The structure of 
the formula for hitting time given by Eq. (3.15) remains the same i.e. 

r/. = r-(Yz,(I-Nz,)-VS)- (3.24) 
3.4.2 Dephasing in the position and coin space 

Complete dephasing has been considered as a kind of decoherence for the discrete-time 
quantum walk on the line [41, 14, 15] and the hypercube in [41]. Here we analyze the 
effect of this kind of dephasing on the hitting time. Consider the discrete-time quantum 
walk on the hypercube. The evolution operator is given hy U = S{I (g) C) where S is 
given by Eq. (2.12) and C is the Grover coin given by Eq. (2.4). The decoherence we 
consider is the dephasing map given by 

V{p) = (1 - p)p + ^pPipPi, (3.25) 

i 

where is a projector onto a basis element of the form \v) \c), v corresponds to a 
vertex and c to a direction in the coin space. The initial condition chosen is |00 . . . 0) (g) 
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{1/y/n) ie., the vertex and the equal superposition of all directions. This initial 

condition lies in the subspace that gives fast hitting times for the non-decohering walk. 
We analyze it in the presence of decoherence. Fig. (3.4) plots the hitting time against the 
decoherence parameter p for dimensions n = 3 to n = 7. Observe that the hitting time 
has an initial jump from j3 = 0top = 0.1 after which it increases smoothly. This jump 
increases from n = 3 to n = 7 suggesting that for higher dimensions it gets worse. Fig. 
(3.5) plots the hitting time for p ranging from 0.01 to 0.1 and we can see that even for 
a little decoherence i.e. p = 0.01, the hitting time jumps and then increases uniformly. 
Fig. (3.6) plots the hitting time against the decoherence parameter for the initial state 
given by |00 . . . 0) (g) |1). We shall show later that this initial state has an infinite hitting 
time for the quantum walk with the Grover coin. As can be seen it becomes finite in 
the presence of decoherence. This shows how this kind of decoherence can turn purely 
quantum behavior into something classical. 

3.4.3 Dephasing in the coin space 

Now we consider a decoherence process that acts only in the coin space. In particular, 
we consider a dephasing map given by 

V{p) = (1 - p)p + Y,p{i ® h)p{i ® h), (3.26) 

i 

where are projectors only onto the coin basis states. Fig. (3.7) plots the hitting time 
against the parameter p. The behavior for this type of decoherence is qualitatively the 
same as before. The hitting time is a little less than in the previous case except at p = 
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Figure 3.4: Hitting time (log scale) vs p for dimensions n = 3 to n = 7 with dephasing in 
the position and coin space. 

and p = 1. At these two points the hitting time is equal to the one in the previous kind 
of dephasing precess. 

3.4.4 Dephasing in the position space 

Finally, we consider a decoherence process that acts only in the position space. The 
dephasing map in this case is given by 



where Pi are projectors only onto the position basis states. Fig. (3.8) shows the variation 
of hitting time with the decoherence parameter. Note that when there is no dephasing 




(3.27) 
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Figure 3.5: Hitting time (log scale) vs p for dimensions n = 3 for small values of deco- 
herence. 

and pure dephasing, i.e. p = and p = 1, then the hitting time is the same whether the 
dephasing is in the position space, coin space or both spaces. 

3.4.5 Analysis of the slope of hitting time 

Here we analyze the slope of the expression for hitting time with respect to the decoher- 
ence parameter p in order to determine the behavior of the jump for small values of p. 
The hitting time from Eq. (3.15) is 




(3.28) 
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Figure 3.6: Hitting time (log scale) vs p for dimensions n = 3 for small values of deco- 
herence for an initial state which has an infinite hitting time. 



where the matrices Y and N are defined as 

Y = {Pf0p}){{i-p)i<S)i + pY^Pi®P*){u<S)U*) 

i 

N = {Qf<E)Q}){{l-p)i<S)i + pY^Pi0p*){U<S)U*). 

i 

Here we consider complete dephasing i.e., Pi is the projector from dephasing in the 
position and coin spaces. Differentiating this expression w.r.t p, we get 
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Figure 3.7: Hitting time (log scale) vs p for dimensions n = 3 to n = 7 with dephasing in 
the coin space. 



The derivative of Y is given by 



^ = (^/^^/)(^)(^^f>*)' (3-29) 



where 

^ = -fc^i + yp.(^P*. (3.30) 

dp 

i 

We can define a similar expression for the derivative of N. In Fig. (3.9) we plot this slope 
for various values of the parameter p. Note that for p = this expression has the same 
problem i.e., of non-invertible matrices as the expression for hitting time. 
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Figure 3.8: Hitting time (log scale) vs p for dimensions n = 3 to n = 7 with dephasing in 
the position space. 

3.5 Prior work 

Kendon [44] has written a good review of the work done so far in analyzing decoherence 
in quantum walks on many graphs such as the line, cycle and the hypercube. Here we 
review some of the prior work on the hypercube. Decoherence in discrete-time quantum 
walks on the hypercube has been studied by Kendon and Tregenna [41], where they 
determine its effect on one-shot and concurrent hitting times defined in [38]. It has been 
proved that the concurrent hitting time is Tc{po) = for pQ = ^( ^^ggi „ ) & hypercube 
of dimension n (for the symmetric initial condition used in this paper). Kendon and 
Tregenna observe that for the range < p < 1/n, where n is the dimension of the 
hypercube, the quantum speed-up is preserved. However, we notice in our simulations 
that for the definition of hitting time considered here (defined in [47]), the situation is 
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Figure 3.9: Slope of the hitting time vs p for the dimension n = 3 with dephasing in the 
position and coin space. 

different. The quantum speed-up is not preserved even for smah values of p, which seems 
to show that this definition of the hitting time is more sensitive to decoherence as it 
captures the time-averaged dynamics. 

Decoherence in continuous-time quantum walks on the hypercube has been analyzed 
in [3] . The kind of decoherence considered is dephasing in the position basis (recall that 
there is no coin space in this kind of walk). The evolution of this walk can be described 
by the Hamiltonian 




where X is the Pauh cj^ operator. The hopping rate is chosen as 7 = k/n, where n 
is the dimension of the hypercube and k is the total energy of the system. In order 
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to incorporate decoherence, a superoperator picture was used. In this picture, the non- 
decohering walk has the superoperator Ut ul, where JJt = exp{—itH). In the presence 
of decoherence this superoperator becomes 

n 

S't = (g) U] (g) • • • (g) (g) • • • (g) [a (g) a], (3.32) 

where 

A = 7i[(lL (g) iX) - {iX (g) 1) - p{i (g) U) + p{Ui (g) Hi) + p{Uo (g) Ho)], (3.33) 

Ho = |0)(0| and Hi = |1)(1|. It has been observed in [3] that for p < 4k the walk has 
quantum behavior and for p > Ak it has classical behavior. 

3.6 Results for the DFT coin 

The hitting time for a discrete-time quantum walk on the hypercube using the DFT 
coin, by contrast to the Grover coin case considered above, can actually be infinite. For 
n=4 we will demonstrate that for the same initial condition as with the Grover coin, the 
hitting time for a quantum walk using the DFT coin is infinity. This is because there 
exist eigenvalues of the evolution operator U whose eigenvectors have an overlap with the 
initial state, but have no overlap with the final vertex for any state of the coin. In this 
subsection we specialize to the DFT coin, but a more detailed and general explanation of 
infinite hitting times based on the symmetry of the walk is presented in the next chapter. 

46 



Suppose there are states \4>) which have no overlap with the final vertex, {(plPflcp) = 0, 
and which are eigenstates of the evolution operator U: U\(f)) = exp(z^)|</)). If the system 
is in the state clearly there is no probability to ever detect the particle in the final 
vertex. Let 'P be a projector onto all such states \4>). Then VPf = PfV = 0, and 
[V, U] = 0, where the commutator is defined as [A, B] = AB — BA. One can write the 
initial state as a superposition of vectors in this subspace and its orthogonal complement: 

\^)=r\^) + ii-r)\^). (3.34) 

Any state that begins in the subspace selected by P will remain there for all time, and 
any state in the orthogonal complement will stay there; this follows from the fact that 
the projectors commute with both the unitary transformation tj and the measurement 
operator Pf. As one starts the walk, the probability that the particle never reaches the 
final state is (*|C/t*P?7*|*), which is 



In order for this probability to be nonzero, there must be eigenstates of the unitary 
evolution operator U which have no amplitude for the final vertex. We can readily 
demonstrate this for the hypercube with the DFT coin. Consider the 4-dimensional 
hypercube. Numerically diagonalizing the evolution operator U, we find it has i, —i, 1 and 
— 1 among its eigenvalues, each with a degeneracy of 8. Since the subspace corresponding 
to the final vertex is 4-dimensional, it is clearly possible to construct a superposition of 
eigenvectors of any of these eigenvalues so that it has no overlap with the final vertex 
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in any coin state. For each of the four degenerate eigenvalues we can construct a 4- 
dimensional subspace of eigenvectors with no overlap with the final vertex, giving a 16- 
dimensional space for all such eigenvectors. By numerically constructing an orthonormal 
basis for this space, we can find an expression for the projector V and measure its overlap 
with the initial state. 

We considered in particular the initial state where the particle was located at the 
|00 . . . 0) vertex and the coin is in an equal superposition of basis states |0), . . . , |n — 1). 
For the hypercube with n = 4 and the given initial state, the probability is 
0.4286, which exactly matches the total probability to never hit the final node after a 
large number of iterations in our numerical simulations. Thus, the probability is close to 
half that the particle never reaches the final state and the hitting time becomes infinity. 

This demonstrates a property of quantum walks not seen in their classical counter- 
parts: for certain initial conditions, there is a nonzero probability that the particle never 
reaches the final state, even though the initial and final states of the graph are connected. 
For a quantum walk with substantial degeneracy, this phenomenon is likely to be generic. 
It might be possible to make the hitting time finite by choosing an appropriate initial 
condition — clearly this happens for the Grover coin — but for some coins this may require 
an initial condition which is not localized on one vertex. Prom our simulations, it it seems 
that for higher dimensions the DFT coin behaves similarly to n = 4. For example, for 
n = 5, our simulations show that the probability to hit the final node increases slowly but 
does not reach 1 even after many time steps. This could be due to the fact that the final 
vertex has no overlap with some eigenvectors of the evolution operator (as for n = 4), 
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and additionally that it overlaps very little for some other eigenvector. This would make 
the probability increase slowly but never reach 1. 

3.7 Results for a quantum walk on a distorted hypercube 

If, as seems likely, the dramatic speed-ups (and slow downs) of quantum walks over their 
classical counterparts depend on the symmetry of the graph, it should be instructive to 
see the effect of deviations from that symmetry. In this section, we look at results for the 
measured walk using the Grover coin on a distorted hypercube. The distorted hypercube 
is defined by constructing the usual hypercube, and then switching two of the connections. 
Pick 4 vertices which form a face-for example, (0 . . . 00), (0 . . . 01), (0 . . . 10), (0 . . . 11). 
Calling these vertices A, B, C, D for short, we distort the hypercube by connecting A to 
D and B to C, and removing the edges between A and B and between C and D. This is 
still a regular graph, and the same quantum walk can be used without having to redefine 
the evolution operator. Unlike the usual hypercube, it is no longer a bipartite graph, and 
the walk can no longer be reduced to a walk in Hamming weight. 

Figure 3.10 plots the hitting time of a quantum walk on a distorted hypercube together 
with that of a classical walk and a quantum walk on regular hypercubes for comparison. 
The hitting time for a quantum walk on a distorted hypercube is more than that of a 
quantum walk on a regular hypercube, but still is much smaller than the hitting time of 
a classical walk. In fact, as the dimension increases one can see that the hitting times of 
the quantum walk on the distorted and regular hypercubes converge towards each other. 
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Figure 3.10: Comparison of hitting times on the regular and distorted hypercubes 

This presumably reflects the fact that for higher dimensions the symmetry is mostly 
unchanged. 



3.8 Infinite hitting times 

We will now show that it is possible for hitting times to be infinite^ and derive a sufficient 
condition for the unitary evolution operator to allow infinite hitting times. As before, 
we begin by forming the projector P onto the subspace spanned by all eigenstates of IJ 
which have no overlap with the final vertex. This projector is orthogonal to the projector 
onto the final vertex, PPj = PfP = 0, and commutes with U, [U, P] = 0. We assume 
(for the moment) that this projector is nonzero; later, we will find a sufficient condition 
for this to be true, and exhibit quantum walks which satisfy this condition. We can write 
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any initial state as a superposition of a state in the subspace projected onto by P and a 
state orthogonal to it, giving the decomposition 



It is easy to see that if |^') Ues entirely inside P, i.e., P\^) = |^), then under the evolution 
the subsequent states will never have any component in the final state and the probability 
defined in eq. (3.6) will be zero. Indeed, since [-P, f/] = and [P,Qf\ = 0, 

p{t) = Tr{PfU[QfUY~^po[mQfY-^mPf} (3.36) 

= T:i{PfU[QfU]^-^PpoP[WQf]^-^WPf} (3.37) 

= Tr{PfPu[QfUf-^po[WQff-^WP^Pf} (3.38) 

= 0, (3.39) 

where po = |^')(^'|. Therefore, the hitting time for this initial state is infinite. More 
generally, if |^) has nonzero overlap with P, P\"^) 7^ 0, then that component of l^*) can 
never reach the final vertex. The probability of ever hitting the final vertex if one starts 
with this initial state is 

p=\{^\{i-Pm\^ <1, (3.40) 

and the hitting time is again infinite. 

To construct this projector, we look at the spectral decomposition of U. If U has 
at least one sufficiently degenerate eigenspace, then we can construct a subspace of this 
eigenspace which has a zero overlap with the final state. For instance, consider one such 
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degenerate eigenspace which has a degeneracy of k. Since the vector space at the final 
vertex is d dimensional (i.e., it has d coin degrees of freedom), we would be solving the 
following d X k system of homogeneous equations: 



V 



N-d+1 



I 



+ 02 



V 



H Vak 



V 



N-d+l 



(3.41) 



Here, we use a labeling where the final vertex in some coin state occupies the last d entries 
of the eigenvectors. The subscript refers to the component of the eigenvector, and the 
superscript distinguishes the eigenvectors in the degenerate eigenspace. This system is 
under-determined if k > d, and it will always have a nontrivial solution — in fact, it will 
have a space of solutions of dimension k — d. Therefore, it is sufficient that there exist 
at least one eigenspace of U with dimension greater than the dimension of the coin, in 
order to have a nonzero projector P. If there is more than one degenerate eigenvalue 
with multiplicity greater than d, the subspace projected onto by P will include all the 
eigenvectors of U which have no overlap with the final vertex. 



The condition derived above is closely related to the question of invertibility of I — N 
in Eq. (3.15) of the previous section. Here we show that I — N is not invertible if and 
only if the projector P is nonzero. Furthermore, in the case when I — N is not invertible, 
the hitting time for a state or density matrix whose support has no overlap with P is 
calculated by replacing the inverse of I — N with its pseudo- inverse in Eq. (3.15). 
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Assume that the projector P is nonzero. Then there is at least one eigenvector \v) of 
U such that Pf\v) = 0. Therefore, 

(Q/f7 ® Q)U*){\v) (8) \v)*) = (Q/ (g) Q)){\v) ® \v)*) = \v) ® 1^;)*, (3.42) 

since tj\v) = exp(i^)|t;) and lJ*\v)* = exp(— and Pf\v) = imphes that Qf\v) = 
\v). Therefore, (I — N)|t;) (g) = 0, I — N has a nonzero nullspace, and hence is not 
invertible. This proves the "if" direction. 

To prove the "only if" direction, assume that I — N is not invertible. This implies 
that there exists a normalized vector \u) eH'SiH* such that 

{i -QfU 0Q}U*)\u) = O^QfU0Qp*\u) = \u). (3.43) 

The vector \u) is an eigenvector of QfU Q*j:U* with eigenvalue 1. Since Qf (g) Q*^ is 
a projector, the vector \u) must therefore lie in the eigenspace of eigenvalue 1 of both 
Qf ® Q*j: and U ®IJ*. This can only be true if \u) is of the form: 

\u) = ^ aij\vi) (g) \vj)*, (3.44) 

where the {l^i)} are eigenvectors of U, and ajj is only nonzero if \vi) and \vj) lie in 
the same eigenspace of U and Pf\vi) = Pf\vj) = 0. Note that the vector \u) need not 
correspond to a physical state in the Hilbert space of the walk. But the existence of such 
a \u) means that the projector P is nonzero, since there must exist at least one 1^,) which 
has a zero overlap with Pf. This proves the "only if" direction. 
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We have shown that if I — N is not invertible then the projector P is nonzero and 
vice-versa. Now we will see that if a density operator p is orthogonal to P, pP = Pp = 0, 
its corresponding vector p^ lies outside the null space of I — N. If p is orthogonal to 
P, then when written in the eigenbasis {\vi)} of U, the diagonal components {vi\p\vk) 
are nonzero if and only if Pf\vi) ^ 0. This implies that for the corresponding vectorized 
quantity p", we have 

r.(P/®P;)pVO. (3.45) 

For any p"", 

r.{Pf ® P}){QfU ® Q/?7*)p" = 0. (3.46) 

Therefore, Np^ / p"" i.e., p"" does not lie in the null space of I — N. Moreover, if p is 
orthogonal to P, then so are J\fp and yM^p. This is easy to see, since from Eq. (3.7) we 
get 

PJ\fp = P{Q fUpU^Qf) = Q fUPpU^Qf = 0, (3.47) 
since PQf = P and [U, P] = 0. Similarly, we obtain Pyp = 0, since 

pyp = P^PfUpU'fPf) = 0. (3.48) 

Therefore, if p is orthogonal to P, then all the terms of the type yN^~^p for all t and 
hence all the terms inside the trace of Eq. (3.9) are orthogonal to P, which means that 
the vectorized versions of all terms inside the trace in Eq. (3.9) lie outside the null space 
of I — N. Thus, for states that do not overlap with P, the hitting time is finite as would 
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be expected, and is given by the same formula Eq. (3.15) with the inverse replaced by a 
pseudo-inverse. 



We focus on the discrete time quantum walk for the remainder of this section, consid- 
ering the walk on the hypercube in particular. It was observed in numerical simulations 
mentioned earlier that for a walk on the hypercube with the DFT coin, an initial state 
given by 



has an infinite hitting time. The phenomenon of infinite hitting times is not restricted 
to the walk with the DFT coin, however. Numerical simulations, followed by analytical 
calculations, have shown that it also occurs with the Grover coin, but for different initial 
states. In fact, it turns out that for the Grover coin, the symmetric initial state (3.49) is 
the only initial state localized at the vertex |00 • • • 0) that has a finite hitting time. Any 
other superposition of the coin states for that vertex will give an infinite hitting time, 
because all such states have a nonzero overlap with P. 

Given any vertex f on a graph, it is natural to ask if there exists any superposition of 
its coin states which overlaps with P, and for which coin state the overlap is maximum 
and for which it is minimum (or zero). We can write the projector P in the form 




(3.49) 




(3.50) 



i,j,k,l 
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where {l^j)} are the vertices and {\k)} are the directions. Suppose the initial state is 



I*) = 1^;) (g)^Q;i|z) = Iv) (g) |a). (3.51) 

i 

Its overlap with the projector P is given by 

mP\^)=^A^^kialai. (3.52) 
k,l 

To find the superposition of coin states such that the overall initial state has the least (or 
greatest) overlap with P, define the matrix 

C"^ = Trvertices |^(|t')(^^| <^ 4om)| , iCv)kl = ^vvkl- (3.53) 

The overlap of the initial state with P can be written in terms of this matrix as 

(*|P|*) = (3.54) 

The matrix Cy is Hermitian and positive, and hence has a spectral decomposition into a 
complete orthonormal basis of eigenstates with non-negative eigenvalues. Assuming that 
{Xi, |ej)} is the spectral decomposition of Cy, we can rewrite the overlap as 

(^|P|^)=^A,|(a|e,)|2. (3.55) 

i 

Prom the above expression, we see that the overlap is maximum (or minimum) if |a) is in 
the direction of the eigenvector with the largest (or smallest) eigenvalue, and zero if \a) is 
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an eigenvector with zero eigenvalue. Therefore, if does not have a zero eigenvalue (i.e., 
is positive definite), then for that vertex every superposition of coin states will overlap 
with P. In other words, the hitting time will be infinity if one starts at that vertex no 
matter what coin state one chooses. Numerical calculations for the Grover and the DFT 
coins on the hypercube show that for the vertex |00 • • • 0), an equal superposition of coin 
states is the only superposition that has a zero overlap with P for the Grover coin, and 
no superposition of coin states has a zero overlap for the DFT coin. Moreover, Tr{P} for 
the Grover coin on the hypercube for n = 4 is 32, which is fully half the dimension of the 
total space (dim=2'^ • 4 = 64). These examples suggest that infinite hitting times may be 
a generic phenomenon on graphs with symmetry. 



3.9 Discussion 

In this chapter, we have examined the definition of hitting time for measured quantum 
walks and analyzed some of its properties. We used this definition to obtain an expression 
for hitting time which is valid on any general graph as long as the unitary evolution 
operator tf of the walk is defined. We simulated this hitting time for a measured quantum 
walk using the Grover coin and compared it to the classical hitting time and to the bounds 
obtained on it; the quantum hitting time is exponentially smaller than the classical hitting 
time. We also showed that the bounds on the hitting time obtained in [38] become less 
tight as the dimension increases. We then investigated the effect of decoherence on the 
hitting time of a discrete quantum walk. We have observed that it affects the hitting time 
more dramatically than it does for the definitions proposed in earlier work. Classical 
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behavior occurs for small values of the decoherence parameter and the hitting time is 
worse than it is for a corresponding classical walk for higher values of this parameter. We 
observed that in the case when the quantum walk has infinite hitting times (a phenomenon 
which is consequence of the quantum behavior of the walk) , decoherence makes the hitting 
time finite and thus make the walk classical. We analyzed the effect of different kinds 
of decoherence and noticed that for dephasing (whether in the position or coin space or 
both), the hitting time deteriorates in more or less the same way. 

These results on fast and infinite hitting times show that simply making a walk quan- 
tum does not guarantee a speed-up over the classical case. We demonstrated that the 
hitting time for quantum walks can depend sensitively on the initial condition, unlike 
classical walks. For certain initial states, the DFT walk can have infinite hitting time, 
a phenomenon not possible in classical random walks. This dependence on the initial 
state varies with the coin used, since for the same initial state the Grover walk has a 
polynomial hitting time. This infinite hitting time is directly related to the degeneracy of 
the eigenvalues of the evolution operator. If the evolution operator is highly degenerate, 
then it is very likely that there exist initial states which give infinite hitting times. 

The cause of the speed-up in quantum hitting time may not be completely clear at the 
moment, but it will be shown in the next two chapters that the symmetry of the graph 
plays a major role in both the speed-up and slow down of the quantum walk. For now, 
we can see this as an interference effect. In the faster quantum walk, the different paths 
leading to the final vertex interfere constructively, enhancing the probability of arrival; 
paths which lead to "wrong" vertices interfere destructively, reducing the probability of 
meandering around in the graph for long times. Unlike a classical random walk, the 
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quantum walk is sensitive to the presence of a global symmetry which is not apparent 
at a purely local level. This phenomenon leads to the speed-up of the continuous-time 
quantum walk on the glued-trees graph as well [19]. 

However, this same reason is undoubtedly the culprit in the slow-down observed for 
the DFT walk. The existence of states which never arrive at the final vertex is made 
possible by the degeneracy of the evolution operator U — a degeneracy which arises due 
to the symmetry of the graph. The existence of states which never arrive at the final 
vertex can also be seen as an interference effect, only in this case the interference of paths 
which lead to the final vertex is destructive: all amplitude to make a transition to the 
final vertex cancels out. 

The connection to symmetry is supported by the quantum walk on the distorted 
hypercube. We observe that the hitting time is worse than that of the usual hypercube, 
but still much smaller than that of a classical walk. The curve of the hitting time on the 
distorted hypercube seems to converge slowly to that of a quantum walk on the regular 
hypercube. This is probably because the distortion used in our simulations is very mild. 
As the dimension grows, and with it the number of edges and vertices, this distortion has 
less effect on the overall symmetry. 

We hasten to add that symmetry of the graph is not the sole reason for speed-ups 
in quantum walks. A polynomial speed-up has been demonstrated in the quantum walk 
versions of the search of an unstructured database [72] and the element distinctness 
problem [4]. However, the dramatic exponential speed-ups have all been demonstrated in 
highly symmetric graphs. In the next two chapters, we explore this idea in more detail 
and show how symmetry can have a strong influence on the quantum walk. 
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Chapter 4 



Symmetry in quantum walks 



We saw in the previous chapter that degeneracy of the evolution operator leads to infinite 
hitting times. One of the main sources of degeneracy in quantum mechanics is symmetry. 
Since the S matrix which makes up part of the evolution operator encodes the connections 
of the graph, it is natural to expect that symmetries of the graph will produce symmetries 
of the evolution operator. Here we analyze this idea. 

Is it sufficient to consider only the symmetries of the graph? Apart from the symme- 
tries induced from the graph, the evolution operator may have additional symmetries of 
its own which lead to additional degeneracy. But such symmetries are difficult to analyze 
in generality, and depend on the details of how one defines the quantum walk. First, 
there is a choice between the discrete and the continuous walk. Second, for the discrete 
walk, if the structure U = S{I (8) C) is used, there is still the freedom to use any unitary 
matrix as the coin. In order to generalize our discussion of symmetries to any walk, we 
restrict attention to the symmetries induced by the graph alone. This naturally leads us 
to the question: "Are there graphs with sufficient symmetry such that, for any walk that 
is defined on the graph, the resultant evolution operator will have enough degeneracy to 
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give rise to infinite hitting times?" It turns out that such graphs do exist, and we give 
an example of such a class of graphs. We will comment briefiy on the effect of additional 
symmetries of the coin for the hypercube later in this chapter. We now look at the au- 
tomorphism groups of Cayley graphs which provide the examples of graphs with infinite 
hitting times. 

4.1 Automorphism groups of Cayley graphs 

An automorphism of a graph is a permutation of its vertices such that it leaves the graph 
unchanged. The set of all such permutations is the automorphism group of the graph. 
When the edge labels or colors in the graph are important, as in the case of a discrete 
quantum walk, we restrict ourselves to those permutations which preserve the edge labels. 
In other words, an edge connecting two vertices has the same label before and after the 
permutation. Such automorphisms are called direction-preserving. In general, we could 
consider automorphisms where we permute the direction labels along with the vertices 
to obtain the same graph with the same coloring. This would form a larger group G of 
which the direction-preserving automorphisms are a subgroup H. 

Since the vertex Hilbert space T-C" has its basis elements in one-to-one correspondence 
with the vertices of the graph, and the coin Hilbert space has a basis in correspondence 
with the direction labels, the automorphisms (which are just permutations of vertices 
and directions) are permutation matrices. In fact, these are all the permutation matrices 
on W (8) H'^ that leave S unchanged, i.e., {all P \ PSP^ = S, where P is a permutation 
matrix}. In this representation, any direction-preserving automorphism has the structure 
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Py®ic, where Py acts solely on T-C" and Ic on H^. Such automorphisms become important 
if we wish to consider the symmetries oi U = S{I (g) C). Clearly, any automorphism of 
this type is a symmetry of C/, since 



(A ® ic) \s{i ® c)] (p„ ® 4)^ = \{Pv ® ic)s{Pv ® 4)^1 {i®c) = s{i ® c). (4.1) 



Elements of G in general do not act trivially on the coin space. Because of this, they 
need not be symmetries of U unless the coin flip operator C respects these symmetries. 



To illustrate all this, consider the example of a hypercube in 2 dimensions (i.e., a 
square). The vertex labels are {(00), (01), (10), (11)} (which also form a basis for T-L"); 
the edges connecting (00) to (01) and (10) to (11) are both labeled 1, and the edges 
connecting (00) to (10) and (01) to (11) are both labeled 2. Thus, the transformation 
(00) ^ (01) and (10) ^ (11), or the transformation (00) ^ (10) and (01) ^ (11), or both 
together, are automorphisms of this graph which need no permutation of the directions. 
Together with the identity automorphism (which permutes nothing) , these permutations 
form the direction-preserving subgroup H. In a matrix representation on the Hilbert 
space TC" ® TC^, they are, 



10 
10 
1 
10 



10 
1 
10 
10 



1 
10 
10 
10 



10 
10 
10 

1 
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where Ic is the 2x2 identity matrix acting on the coin space. These permutations can 
be easily seen to be = {/ O X (g) /, X O / (g) /, X O X (g) /, / O / (g) /}. Just as in the 
representation of S matrix in terms of the Pauh X operators given by Eq. (2.12), this 
group denotes a bit flip in the first, second or both bits of each vertex, together with the 
identity, which gives no flip. (See Fig. 4.1.) 

10 i_ 

2 

00 1 

(00) -«-^(10) 

(01) ^^(11) 

▼ 

00 1_ 

2 

10 i" 

Figure 4.1: The direction-preserving automorphism group of the n=2 hypercube. 

The permutation (10) (01), reflecting along the diagonal while keeping (00) and 
(11) fixed, will be an automorphism only if we interchange the directions \ ^ 2. Sim- 
ilarly, the permutations (00) ^ (11), (00) (01) (11) (10) and (00) (10) 
(11) — > (01) are automorphisms when we interchange the two directions. If we view these 

63 




10 2 00 



11 2 01 



(00) ►-(01) ►-(11) ►-(10) 

11 2 01 




00 2 10 



^(00) ►-(10) ►-(11) ►-(Ol) 

00 2 10 



10 2 00 



Figure 4.2: Automorphisms which interchange directions for the n=2 hypercube. 

permutations along with those obtained above, we obtain a new group G for which H is 
a subgroup. In a matrix representation, the new automorphisms are, 



'^O 1 0^ 



1 
10 



10 



'^l 0^ 



10 
10 



1 



X,, 



^0 10^ 



10 
1 



10 



X,, 



^0 1^ 



10 
10 



10 



Xc, 



where Xc acts on the coin space and corresponds to an interchange of the two directions. 
(See Fig. 4.2.) These four elements of G need not be symmetries of [/, since the coin 
need not be symmetric under conjugation with X^. However, for the hypercube, if we 
use the Grover diffusion matrix as the coin, then the automorphism group G is indeed 
its group of symmetries, since the Grover coin is symmetric under any permutation of 
its basis elements. The symmetry group of the evolution operator would be H if the 
DFT coin is used, since the DFT does not have permutation symmetry. It is important 
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to note that the symmetry group defined above is not the only thing that infiuences the 
degeneracy of the evolution operator. Degeneracy of the coin flip operator C can also 
induce degeneracy in the evolution operator, and the coin may be degenerate even if it 
does not have permutation symmetry (like the DFT coin) . We discuss this in more detail 
later after describing the relationship between symmetry and degeneracy. 

It can be shown that the direction-preserving automorphism group H for any Cayley 
graph is isomorphic to the group on which the graph is defined. This is because any 
direction-preserving automorphism of a Cayley graph is a left translation by a group 
element, and conversely all left translations are direction-preserving automorphisms. The 
first part of the statement is easy to see. Consider any left translation La : G ^ G 
which has the action La{g) = ag, for all g & G. Now, given vertices g and h in G, 
they are connected by an edge from g to h ii g~^h = s, where s E S. Clearly, after 
the transformation we still have {ag)~^{ah) = g~^h = s and hence this automorphism 
preserves the direction labels. Since the group elements are basis states of the vertex 
Hilbert space, a left translation by a group element corresponds to a permutation matrix 
on this Hilbert space. The fact that every direction-preserving automorphism of a Cayley 
graph is a left translation by a group element becomes important in the discussion of 
regular and irreducible representations later in this chapter. 

Finally, let us explicitly construct the representation of the automorphism group for 
the hypercube, which has H = Z2 and G = H ■ Sn- In terms of the Pauli operators the 
representation of H is {ill ■ ■ ■ i^IcXII ■ ■ ■ I^IcXXI ■ ■ ■ i®ic, . . . ,XXX ■ ■ ■ X^Ic}, 
where the tensor product symbol has been dropped in the vertex space, and Ic is the 
identity operator in the coin space. In fact, the representation of H for any Cayley 
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graph will be of the form P 'Si I, where P is a permutation matrix on the vertex space 
and / is the identity on the coin space. The group G for the hypercube will become 
H ■ Sn = {h ■ 7r\h & H,7r G Sn}, where Sn is the permutation group on n elements which 
is assumed to act on the direction labels. 

We now briefly review linear representations of finite groups, and describe how the 
symmetry group of U affects its degeneracy by determining the dimensions of the irre- 
ducible representations of the group. (For further details, see for example [76] or [71].) 

4.2 Representations of finite groups 

A linear representation of a finite group G on a finite-dimensional vector space V is a 
map a : G ^ GL{V) , such that a{gh) = a{g)a{h). GL(V) is the space of invertible 
linear maps of V onto itself. If V is d-dimensional, and we choose a basis of d vectors 
in y, then a{g) for ^ G G becomes a d x d invertible matrix. The trace of this matrix 
is called the character of the representation. Therefore, characters are maps x : G — > C 
with xig) = Tr(c(5'))- Two representations ai and (72 of the group G on vector spaces 
Vi and V2, respectively, are considered equivalent if there exists an invertible linear map 
T : Vi — V2 such that r o ai{g) = (T2{g) o r, for all g E G. The characters of equivalent 
representations are equal, a fact which follows from the cyclic property of the trace 
operator. 

Assume that the vector space V has an inner product defined on it. Since the group 
G is assumed finite, it can be shown that any representation a is equivalent to a unitary 
representation — there exists a basis for V in which a{g) is a unitary matrix for aR g E G 
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(see [71]). A vector space W is said to be invariant or stable under the action of G if 
X G a{g)x G W for all g E G. If the vector space V has a subspace W which is 

invariant under the action of G, then it can be shown that its orthogonal complement 
W-^ is also invariant under G, and so V can be decomposed 



(4.2) 



This means that the representation cr as a matrix on V can be written in a block diagonal 
form consisting of two blocks as 



a\w 



crl 



(4.3) 



where a\w and cr|^± are the restrictions of a to the subspaces W and W-^ 



The linear map a : G ^ GL(V) is called an irreducible representation (irrep), if it is 
a representation and no non-trivial subspace of V is stable under the action of G. Equiv- 
alently, it is an irreducible representation if it is not a direct sum of two representations. 
Any representation on V can, by an appropriate choice of basis, be written as a repre- 
sentation in block diagonal form, where each block corresponds to an irrep. So V can be 
decomposed in the following way: 



V = Wi®W2®---®Wk, 



(4.4) 
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where each of the Wi is stable under the action of cr(G). This decomposition is unique, 
up to reordering of the spaces and an overall equivalence transformation. By an abuse 
of notation, we will use the same labels {V or Wi) to refer both to the vector space and 
to the group representation on that space. Each component Wj of the decomposition in 
(4.4) is isomorphic to an irrep, and the number of such Wj isomorphic to a given irrep 
does not depend on the details of the decomposition. 

If we define an inner product for characters. 



it can be shown that the number of times an irrep with character Xi occurs in a repre- 
sentation with a character % is given by (x|xj) [71]. 

We now define a regular representation of a group G. Suppose |G| = n, and let V 
be an n-dimensional vector space. Let {e^} define a basis for V, which is labeled by 
the group elements of t G G. A regular representation of G is a map a : G V such 
that cr{g)et = Sgt- It can easily be seen that the action of G on the basis vectors is a 
left translation, and in matrix form the representations of the group elements will be 
permutation matrices. An important property of the regular representation is that its 
decomposition into irreps contains all the irreps of the group in it, and each irrep has a 
multiplicity equal to its dimension [76]. Therefore, we can write 




(4.5) 



V = mWi e 712^2 © • • • © rikWk, 



(4.6) 
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where we use the notation nW to mean n copies of the space W, nW = W ®W ® ■ ■ ■ W. 
The {Wi} are all the inequivalent irreps of the group, and rii is the dimension of Wi. 

It will be useful to note that all Abelian groups have one-dimensional irreps. The 
converse of this statement is also true: all groups which have only one-dimensional irreps 
are Abelian. Finally, suppose ai : G ^ Vi and a2 ■ G ^ V2 are two linear representations 
of G on Vi and V2, then ai (g) (J2 : G — Vi (g) V2 is a representation of G on Vi ^¥2- (The 
tensor product (g) of two irreps and of Vi and V2 , respectively, need not 
be an irrep of Vi V2, however.) 

Now consider the unitary operator ?7 on a finite dimensional vector space V which has 
a group of symmetries G. This means that the matrices a{g) representating the elements 
of 5 G G on F all commute with U : [cr(g), U] = 0. Since U is unitary, we can decompose 
V into a direct sum of eigenspaces of U: 

V = Ui(BU2(B---(BUm. (4.7) 
We can also decompose V into a direct sum of irreps of G: 

V = Wi®W2®---®Wk. (4.8) 

It can be shown using Schur's lemma (see [76]) that since G is the group of symmetries 
of U, each irrep of G must lie entirely inside some eigenspace of U. Therefore, for some 
i and j, if Wi C Uj then the degeneracy of Uj is at least equal to the dimension of Wj. 
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Using this fact, we show now that if a graph has sufficient symmetry (in a particular 
sense), then it will lead to quantum walks with infinite hitting times. 



4.3 Discrete-time walks on Cayley graphs 



It was observed earlier that the direction-preserving automorphism group of a Cayley 
graph is the group of left translations of the group elements. Since every group element 
corresponds to a vertex of a Cayley graph and every vertex corresponds to a basis element 
on the vertex Hilbert space (7^") of the walk, the automorphism group of a Cayley graph 
is a regular representation of G on . Every direction-preserving automorphism of a 
Cayley graph will induce a representation on this Hilbert space which looks like a{g) (g) Ic, 
where cr{g) is the regular representation of G on Ti^ and Ic is the identity on the coin 
space. (Note that a{g) (g) Ic is not a regular representation of G.) In order to prove that 
this walk has an infinite hitting time for certain initial states, we need to show that U has 
at least one degenerate eigenspace whose dimension is greater than the dimension of the 
coin. Since every irrep lies completely inside an eigenspace, if one of the irreps occurring 
in (7 (8) J has a dimension greater than the dimension of the coin, then we can say that 
the eigenspace containing that irrep has a degeneracy greater than the dimension of the 
coin. We now show that every irrep of G occurs in cr (g) /. We have. 



Tr{a ® /) = dx, 



(4.9) 
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where % = Tr(cr) and d is the dimension of the coin. If Xi is any irreducible character of 
G, then 

{Tria0i)\x^) = dix\x^)■ (4-10) 

Since x is the character of the regular representation of G, {x\Xi) any irreducible 
character Xi- Therefore, all irreps of G occur in cr®/, and if any irrep of G has a dimension 
greater than the dimension of the coin then there is an eigenspace of U whose dimension 
is greater than the dimension of the coin. So for Cayley graphs, a sufficient condition for 
any discrete-time quantum walk to have infinite hitting times for some initial conditions 
is that the group used to define the graph have an irrep with dimension greater than the 
degree of the graph. 

For a regular graph which is not a Cayley graph, this can be modified as follows. A 
discrete time walk defined on a graph will have an infinite hitting time for certain initial 
states if at least one irrep occurring in the induced representation (on the Hilbert space of 
the walk) of the direction-preserving automorphism group of the graph has a dimension 
greater than the degree of the final vertex. This is a somewhat more difficult to evaluate, 
since unlike a Cayley graph, the induced representation of the symmetry group for a 
general graph is not guaranteed to include every irrep. But in principle it is not difficult 
to check. 

For an example of a graph with infinite hitting times, consider the Cayley graph on 
the symmetric group Sn- T{Sn,X), where X is a generating set for Sn- In order to use 
the form U = S{I ^ C), the graph needs to be |X [-colored and so we chose a generating 
set whose elements x are such that = e, where e is the identity element. For Sn, such 
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a generating set is any set of n — 1 transpositions, e.g., {(1, 2), (2, 3), . . . , (n — 1, n). These 
form the basis for the coin space, and so the dimension of the coin is the cardinahty of the 
generating set: |X| = n — 1. Therefore, a symmetric group Sn which has an irrep with 
a dimension greater than n — 1 will have infinite hitting times for some initial conditions 
for any coin matrix. It turns out that for n > 5 any symmetric group Sn possesses this 
property [36], and so the corresponding Cayley graph r{Sn,X) will have infinite hitting 
times. 

All this indicates that it is not so much the size of the symmetry group that matters 
for infinite hitting times, but rather the kind of group, or more precisely the size and 
number of irreps of the group occurring in the induced representation. Consider the 
hypercube, which has the group H = as its symmetry group (if the coin has no 
permutation symmetry like the DFT coin). This group is Abelian and hence has only 
one-dimensional irreps. So one would expect the unitary evolution operator having this 
as its symmetry group to have very little or no degeneracy. But when the DFT coin is 
used the evolution operator has sufficient degeneracy to have infinite hitting times. This 
is because it is the degeneracy of the DFT coin, rather than the symmetry group, that 
makes the evolution operator degenerate. This is supported by numerical evidence which 
shows that when a randomly generated non-degenerate or slightly degenerate unitary 
coin is used instead of the DFT, the unitary evolution operator has a very small or no 
degeneracy. The degeneracy of the evolution operator with an Abelian symmetry group 
seems to come only from the coin. 
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4.4 Continuous-time walks on Cayley graphs 



In the prior discussion of infinite hitting times we have used the definition (3.5) for the 
hitting time, which is only well-defined for discrete time quantum walks. We have not 
described a suitable measurement process to define the hitting time for a continuous time 
walk; nor is it obvious how to do so in the quantum case, where the presence or absence 
of measurements has a profound effect on the dynamics. Any notion of hitting time for 
the continuous case, however, must include a measurement performed on the final vertex 
at some time which will verify if the particle has arrived there or not. This leaves an 
ambiguity in the definition of hitting time for finite hitting times, but the notion of an 
infinite hitting time still has an intuitive definition: a continuous time quantum walk has 
infinite hitting time if, for any set of measurements on the final vertex at any sequence 
of times, there is always a bounded, nonzero probability that the particle will never be 
found at the final vertex. 

Continuous time quantum walks do not have a coin matrix, and their evolution op- 
erator for undirected graphs is U{t) = ex.p{iHt) where H is the adjacency matrix of the 
graph. Since there is no coin, the degree of freedom at any given vertex is one dimensional. 
Any eigenspace of U with a degeneracy greater than one can therefore contribute to the 
projector P having zero overlap with the final vertex and commuting with U. Whenever 
a measurement on the final vertex is performed, the measurement operators will commute 
with P since P\xf){xf\ = 0, where \xf) is the final vertex state. A nonzero P necessarily 
means an infinite hitting time for initial states that overlap with it. Therefore, only if U 
is completely non-degenerate, and none of its eigenvectors have a zero overlap with the 
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final state, will there be finite hitting times for all initial states. This is a lesser degree 
of degeneracy than is needed in the discrete time case, and we therefore expect infinite 
hitting times to be even more common in continuous time walks than discrete time walks. 

Making the connection to symmetry again, for a continuous walk to have infinite 
hitting times, a sufficient condition is that at least one irrep with dimension greater 
than one occurs in the induced representation of the automorphism group of the graph. 
Consider once more the example of Cayley graphs. As discussed in the case of discrete 
time quantum walks above, the induced representation of the automorphism group G on 
the Hilbert space Ti" is the regular representation. Since there is no coin, Ti" is the Hilbert 
space of the walk. All the irreps of the group appear in this representation because it is 
regular. Therefore, if any of the irreps of G has a dimension greater than one, then the 
walk will have infinite hitting times for certain starting states. Since only Abelian groups 
have all their irreps of dimension one, any Cayley graph defined on a non- Abelian group 
will have infinite hitting times for the continuous walk. 



4.5 Discussion 

In this chapter, we have examined the role of symmetry in shaping the properties of the 
quantum walk. The role of symmetry is not restricted to producing infinite hitting times. 
Symmetry can also be related to the exponentially fast hitting times observed in [38, 47]. 
For one of the cases examined in [47] — the discrete walk on the hypercube with the Grover 
coin — we observe an exponentially smaller hitting time than the classical walk on the same 
graph. But this happens only for the symmetric initial state l^*) = |00 • • • 0) 
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Other superpositions of coin states do not have this speed-up, but rather lead to infinite 
hitting times, because their overlap with P is nonzero. 

We noted that the group of symmetries of the n-dimensional hypercube is G = H ■ Sn 
when one takes into account the direction labels, where H is the normal subgroup of 
direction-preserving automorphisms, and Sn is the permutation group on n elements (in 
this case, the different graph directions), and is also a subgroup of G. It was observed 
in [72] that l^*) is the simultaneous eigenstate of eigenvalue 1 of the subgroup >S„ (more 
precisely, the simultaneous eigenstate of the representation operators of the subgroup). 
Since every element of Sn commutes with U for the walk with the Grover coin, a state 
that begins in an eigenspace of the permutation group will remain in the same eigenspace 
at all times. That is, 

C/*|*) = U'aigm = a{g)U'\^), (4.11) 

where a is the representation of G on the Hilbert space of U and g G Sn- Thus, C/*|^') is 
an eigenstate of a{g) with eigenvalue 1. This eigenspace has dimension 2n. It turns out 
that the final vertex with an equal superposition of coin states |11 • • • 1) (g) Y^- \ i) also 
lies in this eigenspace. Since the walk never escapes this subspace to explore other parts 
of Hilbert space, it leads to an exponentially fast hitting time. For a measured walk, if 
the symmetry subgroup commutes with the measurement operators (as is true in this 
case), then the same argument holds. This shows that those symmetries of the graph 
which are passed on to the evolution operator can create subspaces to which the walk 
may be confined. If the final vertex has no overlap with such a subspace for any coin 
state, then a walk starting in that subspace will have an infinite hitting time. Otherwise, 
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the hitting time will be finite, or even exponentially small depending on the dimension of 
this subspace relative to the full Hilbert space. 

We should also point out that the conditions derived above for infinite hitting times are 
sufficient for a particular graph to have infinite hitting times; but they are not necessary 
conditions. For example, the symmetry group of the hypercube is Abelian, and hence 
does not imply that the evolution operator U must be degenerate. Nevertheless, infinite 
hitting times are observed for quantum walks on the hypercube, due to the fact that the 
choices of coin flip operator (the Grover coin or the DFT) both have their own symmetries, 
which increases the total degeneracy of the evolution operator. Infinite hitting times are 
therefore likely to be even more common than the conditions derived here would suggest. 

One can make very plausible intuitive arguments that both infinite hitting times 
and exponentially fast hitting times are related to symmetry. This makes it seem likely 
that the ideal problem to be solved by a quantum walk would be a problem with global 
symmetry, but in which this symmetry is not apparent at the local scale. In the next 
chapter a framework which explains fast hitting times observed in quantum walks, is 
presented. Using symmetry of the graph again, we show that the quantum walk that 
possesses these symmetries can be confined to a smaller graph called a quotient graph. If 
this new graph is exponentially smaller than the original graph, then it can lead to fast 
hitting times. 
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Chapter 5 



Quotient graphs 



The idea of restricting a search to an invariant subspace of the full search space has proved 
very fruitful in Grover's search algorithm. In both the quantum walk-based algorithm on 
the hypercube in [72] and the "glued-trees" graph in [19], the quantum algorithm works 
very fast by searching a smaller space, where it is known that the solution lies in this 
space. In this chapter, we explore this concept for quantum walks on more general graphs. 
Using symmetry arguments, we show that it is possible to find invariant subspaces of the 
total Hilbert space on which the walk is defined. The automorphism group of the graph 
produces a group of symmetries of the evolution operator for the walk. This group of 
symmetries in turn determines the invariant subspace of the walk. If the initial state is 
in this subspace, the quantum walk effectively evolves on a different graph — a quotient 
graph, which can in some cases be much smaller than the original graph. Here, we give 
a general construction of quotient graphs, given the original graph and a subgroup of its 
automorphism group. We determine the structure of the quantum walk on the quotient 
graph. We then apply the analysis of hitting times developed in [47] and [48] to quotient 
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graphs, and investigate the possibility of both infinite hitting times and reduced hitting 
times on quotient graphs. 



5.1 Action of an automorphism group 

Consider any undirected graph T and let V{T) and E{r) denote its vertex and edge sets. 
Let the graph be colored, not necessarily consistently i.e., the edge between vertices Vi 
and Vj may be colored with a color in the direction from Vi Vj and a color q from 
Vj — Vi- This creates the Hilbert space of positions and colors (or directions) and the 
total space H is spanned by basis vectors \vi,ci), . . . \vn, Cm)- Let this set of basis vectors 
be X. The set of colors at each vertex is not the same for all the vertices since the graph 
may be irregular. We assume that vertices having the same degree have the same set of 
colors. Denote by C„ the set of colors used to color edges going from the vertex v. Thus, 
the shift matrix for this graph is, 



This matrix encodes the structure of the graph F which includes edge colors. An auto- 
morphism of the graph F, as defined above is a permutation matrix which preserves S 
under conjugation i.e., a matrix P such that PSP^ = S. The set of automorphisms form 
a group which we denote by G. 

Now consider a subgroup (not necessarily proper) H of this automorphism group. We 
would like to know what kind of action this subgroup has on the graph and hence on the 
Hilbert space. First, we define what is meant by the term action [67]. 



S = 





(5.1) 
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Definition 5.1. If X is a set and G is a group, then X is a G-set if there is a function 
a : G X X ^ X (called a left action), denoted by a : {g,x) — ^ gx, such that : 

• Ix = X, for all X G X; and 

• g{hx) = {gh)x, for all g,h G G and x E X . 

Definition 5.2. If X is a G-set and x & X, then the G-orbit (or just orbit) of x is 

0{x) = {gx : g e G} C X. (5.2) 

The set of orbits of a G-set X form a partition and the orbits correspond to the 
equivalence classes under the equivalence relation x = y defined hy y = gx for some 
g E G. We can define the action of the subgroup H of the permutation group on the set of 
basis elements X of the Hilbert space H as the multiplication of its matrix representation 
cr{II) (in the basis given by the vectors X) with a basis vector. This is a well-defined 
action since a{l)\x) = \x) and a{g){a{h)\x)) = {a{g)a{h))\x) = a{gh)\x). Therefore, the 
set X is partitioned into orbits under the action of H. 

Since H is a subgroup of the automorphism group, these orbits can be related to the 
graph r through the following results. 

Theorem 5.3. // \v, Ci) and \v, Cj) are in different orbits, then the set of all the vertices 
in the orbits of\v,Ci) and \v,Cj) are the same. 

Proof. If the graph F is irregular (regular graphs are just a special case), then clearly 
any automorphism takes a given vertex to another vertex of the same degree. Thus, 
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automorphisms permute vertices of a certain degree among themselves. Therefore, on 
the Hilbert space H, the matrix representation of any automorphism can be written as 



P = ®Pd, (5.3) 
deD 

where the set D contains all the different degrees in the graph. Consider the subspace 
of vertices of a given degree d which can be written as <8) . Now, if any given 
permutation takes \vi,Ci) to \v2,Cj), then it takes all the basis vectors associated with vi 
to those of V2- Thus, the set of all vertices that lie in the orbit of \vi,Ci) must be the 
same as the set of vertices that lie in the orbit of \vi,Cj) (if \vi,Ci) and 1^1,0^) lie in the 
same orbit, then this is trivially true). Since vi is arbitrary, the set of vertices in the two 
orbits must be the same. □ 

By an abuse of language, say that a vector 1^1,01) is "connected" to \v2,C2) if the 
edge colored ci from vertex vi on the graph is connected to vertex V2 along the color C2 
(i.e., the term \v2,C2){vi,ci\ occurs in S). 

Theorem 5.4. If\vi,ci) and \v2,C2) are "connected," 1^1,01) lies in the orbit Oi and 

\v2-,C2) lies in orbit O2 (not necessarily distinct from Oi), then each of the remaining 
vectors of Oi are "connected" to some vector of 02- 

Proof. If |fi, ci) and |f2, C2) are connected, then there is a term of the type |w2, C2)(fi, ci| 
in S. When we conjugate by some automorphism h E H i.e., perform a{h)Sa{h)^ , then 
this term transforms to c7(/i)|f2, C2)(fi, ci|(t(/i)^ and this must be a term in S because 
a{h)Sa{h)^ = S. This means that the vector that \vi,ci) gets taken to, is connected to 
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the vector that \v2,C2) gets taken to by a{h). Since this is true for all h E H, all the 
vectors in the orbit of 1^1,01) are "connected" to some term in the orbit of |t'2,C2)- D 

The above result applies equally well to any vector in 02- Therefore, one can think 
of the orbits Oi and O2 as being "connected" . 

5.2 Quotient graphs and quantum walks 

Based on the action on a graph F of the subgroup H of its automorphism group, consider 
the following construction of a graph — a quotient graph. The set of vertices occuring in 
an orbit O is a single vertex vq on the new graph and the number of orbits that have the 
same set of vertices is the degree of this new vertex. Thus, for a given vertex vo, the set 
of directions are the various orbits which correspond to the same vertex set. If an orbit 
Oi is "connected" to O2, then the vertices vo^ and vo^ are connected in the quotient 
graph. If Oi and O2 are identical, this corresponds to a self loop for . This means 
that there can be self loops in the quotient graph even if there are none in the original 
graph. We denote the quotient graph obtained by the action of the subgroup on F as 
r/H or Fff. 

Now consider a basis vector \x) = \v,c) and its i7-orbit Ox = {a{h)\x) : h G H}. The 
vector |x) = ^^^^^-^ X^/iei? '''(^)l^) eigenvector of eigenvalue 1 of all the matrices cr(/i) 
forheH since 



aih)\x) = 



1 




1 



^1 h'eH 




(5.4) 
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Similarly, the vector \y) formed from a vector \y) of another orbit is also an eigenvector 
of eigenvalue 1. Each of these vectors {\x)} are orthonormal, since they are formed from 
orbits and distinct orbits do not intersect and they span the simulataneous eigenspace of 
eigenvalue 1 of the matrices cr(iJ). We denote the Hilbert space spanned by these vectors 
by H/H or Hh- Note that the vectors |a;) are just representatives, and any vector in 
its orbit could be used to generate \x). Since the {\x)} are in one to one correspondence 
with the orbits, we let \Ox) denote a vector in H and \x) denote the corresponding basis 
vector in Hh- 

Each basis vector in this space corresponds to a vertex and direction on the quotient 
graph, just as the basis vectors of H, namely {|u,c)} represent a vertex and direction on 
r. Suppose that a vertex v on Fh comes from the set of vertices in the orbit Oi and 
that the orbits O2, ■ ■ ■ ,Ok are all the orbits with the same set of vertices. Since each of 
these orbits is "connected" to some other orbit (either in this set or outside), the degree 
of V is k. Therefore, all the basis vectors lOi), • • • , \Ok) can be associated with v and 
the edges along which they are "connected" to other orbits, as the different directions. 
An alternate labelling of these vectors could be \v,ci), . . . , \v,Ck) and likewise for each 
vertex. Note that this does not produce any natural coloring scheme induced from F, on 
the edges of Th- 

We now show that any discrete quantum walk on F induces a discrete quantum walk 
on Fh as long as U respects H i.e., CF{h)lJa{h)'^ = U \fh e H. Let us define a discrete 
quantum walk as the application of any unitary U which takes a particle on a given vertex 
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V to some superposition of vertices that it is connected to and the directions of only those 
edges which connect them to v. On a basis state it acts as 



U\v,Ci) =^aj\v{cj),c!j), (5.5) 
j 

where |w(cj),c^) and \v,Cj) are "connected", Ylj kiP = 1' ^3 ^'^d the sum runs over all 
the colors of the edges on the side of v. Given this definition for a walk, we have the 
following results. 

Theorem 5.5. Let H he a subgroup of the automorphism group of V and let U be a 

discrete quantum walk defined on T such that U respects the symmetries of the subgroup 
i.e., [U,a{h)] = 0, V/i G H. If the initial state lies in the subspace spanned by all the 
orbits {\Oi)} under the action of H, then the walk is contained in the subspace. 

Proof. We have 

U'\Oi) = U'a{h)\Oi) = a{h)U'\Oi). (5.6) 

This shows that since \Oi) lies in the eigenspace of eigenvalue 1, U^\Oi) also lies in the 
same space, which is spanned by {|Ci)}. □ 

Theorem 5.6. Let H be a subgroup of the automorphism group of T and let U be a 

discrete quantum walk defined on F such that U respects the symmetries of the subgroup 
i.e., [tJ,a{h)] = 0, G H. If the initial state lies in the subspace spanned by all the 
H-orbits {\Oi)}, then U induces a walk on Fh in the Hilbert space Tin- 
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Proof. In order to show that U induces a walk on Th, we need to show that its action is 
similar to Eq. (5.5): 

U\0) = Y,bj\Oj), (5.7) 
j 

where \Oj) are the orbits connected to \0). But this follows from the fact that if the walk 
moves the particle from a vector to vectors "connected" to it, then it does the same for 
superpostions of vectors or the orbit states \0). □ 



We can derive the structure of this induced walk from the original walk by making 
use of its action on the orbit states. The induced walk on the subspace Hh becomes 
Uh = y{^y\^\^x)\y){^\- This defines a unitary operator in Hh because, 

= Yl{^y'\U^PHU\Oy)\y'){y\ 

i,y,y' 

= Ih, (5.8) 

since {/ commutes with Ph, where Ph is the projector onto Hh- Now consider the shift 
matrix of the walk. Its action on Hh is given by Sh = J2x yi^y\'^\^x)\y)i^\- '^^^ 
expression {Oy\S\Ox) is non-zero if and only if the two orbits are "connected". If two 
orbits are connected then they must be a superposition of the same number of vectors 
i.e., \Ox\ = \Oy\ and each vector in the superposition in \Ox) is connected to one vector 
in the superposition in \Oy). Therefore, 



{Oy\S\Ox) = \Ox\/^\Ox\\Oy 
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Thus, 

SH = Y.\y){x\. (5.9) 

x,y 

This means that the action of Sh is very similar to the action of S in that it takes the 
walker from any vertex to the vertex it is connected to in the quotient graph. The action 
of the coin which was / (g) (7 on the original graph becomes Ch on the quotient graph so 
that Uh = ShCh- Moreover, Ch can be decomposed as follows, 

Ch = Ci®C2®---®Cn, (5.10) 

where N is the total number of vertices of the quotient graph and each Q acts only on 
the basis vectors associated with the vertex Vi of the quotient graph and each Ci has a 
dimension di which corresponds to the degree of the Vi. In the following examples, such 
a decomposition is provided along with a list of the basis vectors on the quotient graph 
such that Ci acts on the first di basis vectors, C2 acts on the next ^2 vectors etc. 

5.3 Examples of quotient graphs 

In this section, we illustrate the above abstract formalism with some examples. In all 
of the examples we use the following notation to describe the subgroups used to find 
quotient graphs. The elements of the subgroups denote permutations of directions, but 
it is to be understood that this has to be done along with an appropriate permutation 
of vertices, which makes it an automorphism of the graph. Although such a permutation 
of vertices need not exist for every permutation of directions, they exist for the examples 
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that we consider here. Moreover, this permutation of vertices can be specified simply: 
permute the generators which are in one-to-one correspondence with the directions in the 
same way as the directions and this induces a permutation of vertices. 

For example, let (1,2) be a group element. This is the automorphism obtained by 
interchanging directions 1 and 2 and interchanging generators ti and t2 so that vertices 
such as tit2 go to t2ti etc. We do not consider direction preserving automorphisms in 
the following examples, since they tend to give rise to quotient graphs with self loops. 
Finally, we use cycle notation to denote permutations, i.e., (1,2,3) means 1 goes to 2, 2 
goes to 3 and 3 goes to 1. 

Example 1. As the first example, consider the Cayley graph r(S'3, {(1, 2), (2,3)}), and 
let {ti,t2} = {(1, 2), (2, 3)}. The basis vectors of the Hilbert space of the walk are 
{|e, 1), |e, 2), 1), . . . , |tii2^i, 2)}. The automorphism group of this graph is Aut(r(6'3, T)) ~ 
R{S3)Z2. Consider the subgroup H = Z2 which corresponds to interchanging the direc- 
tions 1 and 2. The orbits under the action of this subgroup are 



(l/V2)(|e,l) + |e,2)), 



O2) 



(1/V2)(|ii,l) + 1*2,2)), 



O3) 



(l/V2)(|ii,2) + |i2,l)) 



(1/V2)(|M2,2) + 1*2*1,1)), 



(1/V2)(M2,1) + Mi,2)), 



(l/V2)(|ii*2*i,l) + |M2ii,2)). 
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The original and the quotient graph in this case are shown in Fig. 5.1. The unitary 
describing the quantum walk on F is given hy U = 5(7 (g) (7) where S = S' + S'^ and 

S' = \e, l){h, 1| + \e,2){t2, 2| + {h, 2){ht2,2\ + \t2, l){t2ti, 1| + 1*1*2, l){tit2h, 1| 
+ Mi,2)(Wi,2|. 

This becomes Sh on the quotient graph and is given by Sh = S'h + and 

S'h = \0i){02\ + \03){04\ + \05){0e\. (5.11) 

This can also be written by giving new labels to the vertices and directions of the quotient 
graph: 

S'h = \vi,R){v2,L\ + \V2,R){V3,L\ + \v3,R){v^,L\, (5.12) 

where we have relabeled \Oi) through \Oq) as \vi,R) through \v4,L). Note that there is 
no \vi,L) and \v4,R) which exactly corresponds to the way these vertices are connected 
in the quotient graph. Now, if we take the coin to be C = cTa;, the Pauli X operator 
(which is also the Grover coin in two dimensions), then on the quotient graph the coin 
flip matrix Fh = {I ®C)h becomes 

Fh = \vuR){vuR\ + \V4, L){v4, L\ + F'h + F'^, (5.13) 
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where = \v2, L){v2, R\ + {v^, L){v3, R\. It can also be written as 



Fh = i(BC' (BC' ei, 



(5.14) 



where C' = X, the Pauh ax operator. Thus, the walk becomes Uh = ShFh i.e., 



Uh 



( \ 

1 

1 
1 
1 
1 



1 





Quotient graph under H 



t^tg 



Figure 5.1: The graph r(S'3, {(1, 2), (2, 3)}) and its quotient graph. 
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Example 2. Now consider the Cayley graph r(53, {(1, 2), (2, 3), (1, 3)}) where {(1, 2), (2, 3), (1, 3)} 
= {ti,t2,t3}. A subgroup of its automorphism group is Ss which consists of all per- 
mutations of the three directions. Consider a subgroup of this consisting of Hi = 
{e, (1, 2, 3), (1, 3, 2)}. Under the action of this subgroup, the orbits are 

= (l/x/3)(|e,l) + |e,2) + |e,3)), 

\02) = (l/\/3)(|ti,l) + |t2,2) + |t3,3)), 

\Os) = (l/V3)(|ti,3) + |i2,l) + K3,2)), 

104) = (l/\/3)(|ti,2) + |i2,3) + 1*3,1)), 

105) = (l/\/3)(|tit2,l) + 1*1*2,2) + |tlt2, 3)), 

\Oe) = (l/V3)(|t2ti,l) + |*2ti,2) + |i2ti,3)). 
The shift matrix for this walk becomes Shi = S'jj^ + , where 

S'h, = \0i){02\ + \0e){03\ + |05)(04|. (5.15) 
We can relabel the quotient graph as shown in Fig (5.2). The matrix S'^j^ becomes 

S'hi = \n, 1){V2, 1| + |^^3, ^){V2, 2| + \V4, 1){V2,3\. (5.16) 
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If we choose the Grover coin for the walk, the walk on the quotient graph becomes 



0—12 2 Q Q 

u 3 3 3 



1 



1 



1 



4 



^ 

" 3 







(5.17) 



Fig. 5.2 also shows the quotient graphs for the above Cayley graph with subgroups H2 ^ 
^3 and Hs ~ {e, (2, 3)}. The basis states of the quotient Hilbert space Hh2 are 



\Oi) 



\02) 



(|e,l) + |e,2) + |e,3))/V3, 



(|ti,l) + |t2,2) + |t3,3))/V3, 



103) = {\h,2) + \h,3) + \t2,l) + \t2,S) 
+ 1^3,1) + 1^3, 3))/V6, 

104) = (|iii2,l) + |iii2,2) + |tit2,3) + |t2ti,l) 
+ \t2ti,2) + \t2ti,3))/V6, 
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and the basis states of Hhz are 



Pi) 


= Kl), 




= (|e,2) + |e,3))/V2, 


Ps) 






= (|ti,2) + |ii,3))/V2, 




= (1^2,2) + |t3,3))/V2, 


\ 

1^6/ 


— I\+ 'i\ _l_ If 9\'\ /a/9 




= (1*2,1) + 1*3, 1))/V^, 


Ps) 


= (1*1*2,3) + 1*2*1, 2))/V2, 


Pq) 


= (1*1*2,1) + 1*2*1, 1))/V2, 


Pio) 


= (1*1*2,2) + 1*2*1, 3))/V2. 



The unitary corresponding to the walk on the quotient graph of H2 



-1/3 2^2/3 0^ 



1 

-2V2/3 1/3 Oj 



and the one on the quotient graph of H3 is 




Figure 5.2: The graph r(53, {(1, 2), (1, 3), (2, 3)}) and its quotient graphs. 



The matrices Sh^ and Chs are given by, 5^3 = 5" + 5""^ and Chs = C" C" © C" © C" 
where, 

5' = {\Oi){03\ + \02)m + \04){Oio\ + \OQ){Og\ 

+ mm, 



c' 



'^-1/3 2y/2/^ 
^2^/2/3 1/3 ^ 
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and 



^ -1/3 2/3 2/3 ^ 



2/3 -1/3 2/3 
2/3 2/3 -l/3y 



Example 3. In this example, we determine the quotient graph of r(S'4,r) for T = 
{(1, 2), (1, 3), (1, 4)} = {ti,t2,i3} under the subgroup H Ss which corresponds to aU 
possible permutations of the directions at each vertex. The original and the quotient 
graphs are shown in Fig (5.3), where "t" has been dropped in the vertex labels. There 
are 14 orbits under the action of this subgroup. These are 

\Oi) = (|e,l) + |e,2) + |e,3))/V3, 
IC2) = (|ti,l) + |t2,2) + K3,3))/V3, 

103) = (|il,2) + |i2,l) + |t3,2) 

+ \t2,S) + \ti,3) + \t3,l))/V6, 

104) = {\tlt2,l) + \t2tl,2) + \t3t2,3) 

+ \t2t3,2) + \t3ti,3) + \tit3,l))/V6, 
IC5) = {\ht2,2) + \t2h,l) + \t3t2,2) 

+ \t2t3,3) + \t3ti,l) + \tlt3,3))/V6, 
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\0&) = (|tlt2,3) + |t2tl,3) + |i3t2,l) 

+ \t2h^) + \hh,2) + \tlh,2))/^/l, 

IO7) = {\tlt2tl,l) + \tlt2h,2) + \t^t2h,2,) 

+ 1*3*2*3, 2) + 1*1*3*1, 3) + 1*1*3*1, 1))/V6, 

Ps) = (1*1*2*1, 3) + 1*3*2*3, 1) + 1*1*3*1, 2))/V^, 

IO9) = (1*1*2*3,1) + 1*2*1*3, 2) + 1*3*2*1, 3) 

+ 1*1*3*2, 1) + 1*2*3*1, 2) + 1*3*1*2, 3))/V6, 

|Clo) = (1*1*2*3, 2) + 1*1*3*2, 3) + 1*2*1*3,1) 

+ 1*2*3*1, 3) + 1*3*1*2, 1) + 1*3*2*1, 2))/V6, 

Pll) = (1*1*2*3,3) + 1*1*3*2,2) + 1*2*1*3,3) 

+ 1*2*3*1, 1) + 1*3*1*2, 2) + 1*3*2*1, 1))/V6, 

IC12) = (1*3*1*2*1, 1) + 1*2*3*2*1, 3) + 1*3*1*2*1, 2) 

+ 1*2*3*2*1, 2) + 1*1*3*1*2, 1) + 1*1*3*1*2, 3))/V6, 

IO13) = (|ti,l) + |i2,2) + |t3,3))/^/3, 

IC14) = (1*1*3*2*1, 1) + 1*1*3*2*1, 2) + 1*1*3*2*1, 3) 

+ 1*2*3*1*2, 1) + 1*2*3*1*2, 2) + 1*2*3*1*2, 3))/V6. 

The unitary walk on the quotient graph can be written as 

Uh = Sh- Ch. (5.20) 
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The matrices Sh and Ch are given by, Sh = S'+S'^ and Ch = 1®C'®C"®C'®C"®C'®1 
where, 

S' = {\Oi){02\ + \03){0^\ + \05){07\ + \Oe){Og\ 

+ \08){Oi2\ + |Ol0)(Ol3| + |Oll)(Ol4|), 



and 



^ -1/3 2^2/3^ 
^2^/2/3 1/3 ^ 



-1/3 2/3 2/3 ^ 
2/3 -1/3 2/3 
\^ 2/3 2/3 -1/3 J 



Example 4. Consider the hypercube. The automorphism group of the hypercube is 
Aut(r(2^2'! ^)) — ^2^n- We focus on the subgroup Hi = Sn and look at the resulting 
quotient graph. We consider the case when n = 3, but the procedure for a general n is 
very similar. The subgroup Hi consists of all possible permutations of n directions. The 
orbits under the action of this subgroup are given by 

\Oi) = (|000,1) + |000,2) + |000,3))/V3, 

102) = (|001,1) + |010,2) + |100,3))/V3, 

103) = (|001)(|2) + |3)) + |010)(|1) + |3)) 
+ |100)(|1) + |2)))/V6, 
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Figure 5.3: The graph r(S'4, {(1, 2), (1, 3), (1, 4)}) and it quotient graph under the sub- 
group H. 



104) = (|011)(|1) + |2)) + |101)(|1) + |3)) 
+ |110)(|2) + |3)))/V6, 

105) = (|011,3) + |101,2) + |110,1))/V3, 
lOe) = (|111,1) + |111,2) + |111,3))/V3. 



The graph becomes a Hne as shown in Fig. 5.4 and all the vertices of a certain Hamming 
weight collapse to a point. This fact h first been observed in [57]. In [72], this idea 
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was used to construct a search algorithm on the hypercube. As observed in [72], the 
states on the Hne can be relabeled |0, R) L) R) , |2, L), |2, R), |3, L). For the general 
hypercube of dimension n, these states generalize to 



\x\=xxa=0 



1^'^) = JyiE (5.21) 



where |x| is the Hamming weight of x. 





000 



100 
010 
001 



110 

101 

oil 



Figure 5.4: The n = 3 hypercube and its quotient graphs. 
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Under the action of a different subgroup H2 = Sn-i consisting of permutations of 
n — 1 directions and the corresponding permutations of vertices, the quotient graph is 
shown in Fig. (5.4). The basis states of 'Hh2 when n = 3, are 





= |000,1), 




= (|000,2) + |000,3))/V2, 


Ps) 


= |ooi,i), 




= (|001,2) + |001,3))/V2, 


\05) 


= (|010,2) + |100,3))/V2, 




= (|010,1) + |100,1))/V2, 




= (|010,3) + |100,2))/V2, 




= (|011,2) + |101,3))/V2, 


P9) 


= (|011,1) + |101,1))/V2, 


Pio) 


= (|011,3) + |101,2))/V2, 




= (|110,2) + |110,3))/V2, 


P12) 


= |iio,i), 


P13) 


= (|111,2) + |111,3))/V2, 


\Ou) 





For any general n, the graph is still planar as shown in Fig. 5.5 and there will be 6n — 4 
basis states. They can be labeled as \xO, L),\xO, R),\xO, D),\xl, L),\xl, R),\xl,U), where 
X is the Hamming weight of the last n — 1 bits (which fall under the action of the subgroup 
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xO,x=0 



xO,x=1 



x0,x=2 



xO,x=n-1 



x1 ,x=0 



x1 ,x=1 



x1 ,x=2 



x1 ,x=n-1 



Figure 5.5: The quotient graph of a general hypercube under the group Sn-i- 



Sn-i) and the bit next to it is the first bit. L,R, U and D stand for left, right, up and 
down respectively. They are given by 



l-l-«> = -.)("-') ,1^ ^J"'-'*' 

'\ X ! \x\=xXci=Q 



\x\=x Xcl= 



\x\=xxa=l 



|x|=x 
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Note that the states \xO,U) and \xl,D) do not exist. Moreover, \xO,L) and \xl,L) do 
not exist when x = and \xO,R) and \xl,R) do not exist when x = n — 1. The unitary 
matrices describing the walk on these graphs are 



and 



-1/3 2^2/3 

1 
-1/3 2^2/3 
2^2/3 1/3 
1 
2V2/3 1/3 



(5.23) 



(5.24) 



The matrices and are given by, = S' + S'^ and = C" C" C" C" 1 
C' C" where. 



S' = {\Oi){03\ + \02){05\ + \Oi){Os\ + \OG){09\ 



+ p7)(Cll| + |Cio)(Ol3| + |Cl2)(Cl4|), 



-1/3 2^2/3^ 
^2^2/3 1/3 ^ 



100 



and 



1/3 2/3 2/3 



C' 



2/3 -1/3 2/3 



2/3 2/3 -l/3y 



Example 5. While we have shown how to construct quotient graphs for discrete-time 
walks on Cayley graphs, the idea of a quotient graph is more general. In this example, 
we consider the "glued trees" graph shown in Fig. 5.6. This graph is not regular and 
hence not a Cayley graph. It is undirected, and so we can easily define a continuous 
walk on it. Because the continuous walk does not have a coin space, we need not consider 
permutations of directions in the automorphisms. Quantum walks on this graph were first 
analyzed in [18], and it has been shown that quantum walks move exponentially faster 
on this graph from "entrance" to "exit" than classical walks. The main reason for this 
exponential speed up is that the quantum walk moves in a superposition of all the vertices 
in a given column. It can be seen that in any given column, the vertices which branch 
out from the same vertex in the previous column can be interchanged as long as the 
corresponding interchange on the other side of the central column takes place. Therefore, 
the automorphism group of this graph is where k is one half of the total number of 
vertices on one side of the central column. Under the action of these automorphisms, the 
vertices in each column form a single orbit, and hence collapse to a single point in the 
quotient graph. There are 2n + 1 orbits under the action of this subgroup, where the 
columns j are such that < j < 2n. The orbits can be written as 





(5.25) 



vG column j 



101 



The Hamiltonian for the quantum walk on the quotient graph becomes [40] 



{j\H\j±l) = -V2j 

27 j = 0, n, 2n 
37 otherwise, 



m\~j) 



(5.26) 



with all other matrix elements zero. This is also shown in Fig. 5.6 where the 7 has been 
dropped for brevity. 




-V2 -V2 -V2 -V2 -V2 -\f2 -V2 -V2 
• • • • • • • • • 

233323332 

Figure 5.6: The glued trees graph and its quotient graph. 
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5.4 Automorphism group of the quotient graph 

In this section we determine the automorphisms of the quotient graph which are induced 
from the automorphism group of the original graph and the subgroup used to obtain 
the quotient graph. In deahng with the automorphisms of T we used permutations of 
vertices and edges, and this in turn corresponds to permutations of basis vectors which 
preserves the shift matrix. On the quotient graph, we define those permutations of orbits 
which preserve the Sh matrix as automorphisms, since there is no natural choice of edge 
colors. These permutations of orbits which preserve the new shift matrix also preserve 
the quotient graph. 

Let Gi be a set of automorphisms of T which are of the following type. If they take a 
basis vector belonging to a iJ-orbit Oi to a basis vector belonging to O2, then they take 
every basis vector in Oi to some basis vector in 02- Clearly, all the automorphisms in H 
are of this type, under the special case when Oi = 02- It is also easy to verify that Gi 
is a subgroup of G, and that H is a subgroup of Gi. We now show that H is a normal 
subgroup of Gi i.e., ghg"^ e H,'ig e Gi and V/i G H. 

Theorem 5.7. Given the group Gi defined as above we have, 

1. The subgroup H is a normal subgroup of Gi. 

2. G\ is the largest subgroup of G such that H is a normal subgroup of Gi — that is, 
for any g E G, if ghg"^ E H E H then g E Gi. 

Proof. We show this by considering the action of all of these group elements on the set 
of basis vectors. 
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1. Let X be any basis element belonging to some iJ-orbit Oi, and let g take every 
element in O2 to some element in Oi. Then ghg~^x = ghy, where y G 02- Now, 
hy = z where z e O2 since these orbits are formed under the action of H. Hence, 
gz = x', where x' e Oi. But every x' G Oi can be written as h'x for some h' G H. 
Thus, ghg-'^ G H. 

2. Consider some basis element x e O3 and let gx G O4. Since ghx = h'gx, g{hx) = 
h'y = y', where y,y' G O4. Therefore, g takes hx G O3 to y' G O4, but since h e H 
is arbitrary, 3 takes every element of O3 to some element of O4. It follows that 

□ 

Since H is normal in Gi, the quotient set Gi/H, i.e., the set of all cosets gH = 
{gh\g G H}, is a group. This group has a natural representation in the Hilbert space Hh 
as a permutation matrix in the basis where each orbit is a basis vector. 

Theorem 5.8. Gi/H c Aut{TH). 

Proof. Consider any automorphism g G Gi and let a{g) be its representation in H. Then, 
we have a{g)Sa{g~^) = S. The projection of this into Hh is given by 



PHaig)Saig-')PH = PrSPr = Sh- 



(5.27) 
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The representation a{g) commutes with Ph, since it permutes all the vectors in an orbit 
to vectors in another orbit. Therefore, 

PHCT{g)PHSPHa{g-^)PH = Sh- (5.28) 

But as a representation, PH(T{g)PH = a{gH). This means that the representation of 
gH e Gi/H in Hh is a group of symmetries of Sh and therefore G\/H (Z Aut(r). □ 

We see that the quotient graph is obtained from F modulo the symmetries in H. 

5.5 Hitting time on quotient graphs 

In this subsection, we address the question of when quantum walks on quotient graphs 
have infinite hitting times. It is possible that for some subgroups, the walk on the quotient 
graph does not have infinite hitting times even if the walk on the original graph does. 
In order to carry over the discussion of hitting times to quotient graphs, we must keep 
in mind that the evolution operator is now followed by a measurement. To remain on 
the quotient graph (i.e., in the subspace given by Ph), the measurement operators must 
commute with the symmetry operators (j{h), h E H. 

If this condition is satisfied, then we can obtain a condition to check whether the 
quotient graph has initial states with infinite hitting times: if the subspace of those 
initial states with infinite hitting times on the original graph whose projector is P, has no 
nontrivial intersection with the subspace whose projector is Ph i.e., PDPh = 0, then the 
walk on the quotient graph does not have infinite hitting times. If there is a nontrivial 
intersection, then it does. (Here we have used the projectors onto the subspaces to denote 
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the spaces themselves.) This condition can also be verified by obtaining the restriction 
of the evolution operator and the measurement operators onto the quotient graph. By 
diagonalizing the new unitary evolution operator and constructing the projector of states 
P which have no overlap with the new final vertex state. The subspace of these states is 
exactly the intersection PHPh- We will examine this condition for some of the examples 
considered above. 

In the first example, we choose the final vertex to be tit2ti. The measurement op- 
erators are Pf = \t1t2t1) {tit2ti\ (8) I and / — Pf. This measurement commutes with the 
subgroup chosen, and the quotient graph does not have infinite hitting times. This is 
because the original graph does not have infinite hitting times either i.e., P = 0. 

For the second example, for the graph r(5'3, {(1, 2), (2, 3), (1, 3)}), we used three dif- 
ferent subgroups and form their quotient graphs. In order to determine whether the 
quotient graph has infinite hitting times for various subgroups, we must choose different 
final vertices for the different subgroups since the measurement must commute with the 
symmetries. Therefore, for Hi we choose tit2 as the final vertex, and the measurement 
operators are 1*1*2) (^1*21 <X) / and its orthogonal complement. This measurement com- 
mutes with the subgroup Hi. For this final vertex and measurement, the original graph 
has infinite hitting times i.e., P 7^ and the quotient graph also has infinite hitting times 
i.e., PCiPh^ / 0. In fact, using the C-matrix defined above in Eq. (3.53), we find that if 
the initial vertex is the identity |e), then there is no superposition of coin states that has 
a finite hitting time, because C„ does not have a zero eigenvalue for v = e. 

For the subgroup H2, we choose the final vertices to be titg and t2*i- Therefore, 
the measurement on the original graph must be a projective measurement with outcomes 
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— (I(iii2)(iii2| + |i2ii)(i2*i|)®-f and its orthogonal complement. For this measurement 
and final vertices, the original graph has P = 0. Therefore, the quotient graph also does 
not have infinite hitting times. For the subgroup H^, the measurement operators are 

— (I(iii2)(iii2| + |*2ii)(i2ii|)<^-f and its orthogonal complement. For this measurement, 
neither the original graph nor the quotient graph have infinite hitting times. 

In example 3, for the Cayley graph r(S'4, {(1, 2), (1, 3), (1, 4)}), we choose the final 
vertices to be \t1t3t2t1) and \t2t3tit2)- In this case, we find that while the original graph 
has infinite hitting times, the quotient graph does not. In fact, on the original graph, 
the equal superposition of all coin states at the vertex |e) is the only superposition which 
does not have an infinite hitting time (i.e., the Cy matrix has only one zero eigenvalue 
with the equal superposition of coin states as its eigenvector). It is precisely this vector 
which is included in the subspace of the quotient graph. This is not a coincidence — in 
both cases, it is picked out by the symmetries of the graph. 

In example 4, for the two different subgroups of the automorphism group considered 
for the hypercube, we find that the behavior of hitting times is very different. For the 
subgroup Hi , the quotient graph becomes a line with the vertex = 00 ... on one end 
and the vertex 1 = 11 ... 1 on the other. If one designates the final vertex to be 1 by 
choosing the measurement operators to be Pf = |1)(1|(8)/ and its orthogonal complement, 
then we find that this quotient graph does not have infinite hitting times for any initial 
state: PDPhi = 0. In fact, if the initial state is |00 . . . 0) (8) ^ J2i K)' ^^^^^ the hitting time 
is polynomial in d, the dimension of the hypercube [38, 47]. Using the C-matrix for the 
original graph, we find that if the initial vertex is 00 ... 0, then the equal superposition 
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of all directions is the only zero eigenvector of Cy, which means that it is the only coin 
state that does not have an infinite hitting time. 

On the other hand, choosing 11 ... 10 as the final vertex and using the subgroup H2, 
we find that the quotient graph shown in Fig. 5.4 does have infinite hitting times for some 
initial states i.e., PflPffj / 0- Using the C-matrix again, we find that if the initial vertex 
is 00 ... 0, then the equal superposition of all directions once again is the only coin state 
that has no infinite hitting times. For every other superposition of coin states for that 
vertex (i.e., every other \a) in Eq. (3.51)) Cy has a nonzero eigenvalue. 



5.6 Discussion 



We have investigated the behavior of quantum walks on undirected graphs by making 
use of the automorphism group of the graph. Automorphisms of the graph may become 
symmetries of the discrete quantum walk, depending on the symmetries of the coin matrix. 
Quantum walks which respect the symmetries of some subgroup H of this automorphism 
group have an invariant subspace in the total Hilbert space. We showed that the walk 
restricted to this subspace can be seen as a (different) quantum walk on a quotient graph, 
and that this graph can be constructed from the original graph given the subgroup H. The 
dynamics of the new walk can also be derived from the original walk and the subgroup. 
The quotient graph is obtained from the original graph by identifying vertices and edges 
which form an orbit under the action of H; this means that the quotient graph and the 
new quantum walk both have no symmetries coming from H. The new quantum walk 
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only has the remaining automorphisms as its possible symmetries, and so it has, in a 
sense, "used up" the ones in H. 

To discuss hitting times, we use the measured walk defined in [38] and [47], which 
consists of the application of a unitary operator followed by a projective measurement at 
each time step. For the walk on the quotient graph to be preserved, the choice of mea- 
surement must commute with the symmetries in H. This restriction is very important: 
even if the walk and initial state both have a larger group of symmetries, the walk will 
be on a quotient graph corresponding to a smaller subgroup H if the measurement does 
not commute with the remaining elements of the larger group. 

For instance, in Example 2, using the subgroup Hi, we obtained a walk on its quotient 
graph. Suppose the measurement is a projective measurement of the vertex tit2- The 
initial state je) (g) (jl) + 12) + 13))/-^ (3), and the walk with the Grover coin U both have all 
the symmetries of H2. But the measurement, which commutes with all the elements of 
Hi, does not commute with all those in H2, and the effective walk will be on the quotient 
graph corresponding to Hi. 

The remaining symmetries of the evolution operator can lead to degeneracy in its 
eigenspectruni [48], and may result in infinite hitting times on the quotient graph. In 
general, we found a condition to determine whether the walk on the quotient graph will 
have infinite hitting times: given the original graph, the quantum walk and any subgroup 
H, one can determine the projector onto states with infinite hitting time P, and the 
invariant subspace of the quotient graph Ph- If PDPh = then the quotient graph does 
not have infinite hitting times. 
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Even when the hitting time is not infinite for an initial state on the quotient graph, it 
is possible that it could be extremely long. It would be useful to have a criterion to pick 
out subgroups of the automorphism group whose quotient graphs have exponentially fast 
hitting times. For example, in the case of the hypercube, using the subgroup Hi (whose 
quotient graph is a line) turns out to give very fast hitting times. But on a general 
undirected graph it is not easy to determine whether there is a subgroup whose quotient 
graph gives fast hitting times. 

To investigate this, we need to make the notion of "fast" more precise. One way to 
define "fast" for a parametrized class of graphs (such as the hypercube, where the pa- 
rameter is the dimension) is to say that the hitting time must be 0(log N) (exponentially 
smaller), where A'' is the number of vertices of the graph. Using this notion, we can expect 
fast hitting times to exist in graphs which have quotient graphs with an exponentially 
smaller number of vertices. While this is not necessarily a sufficient condition for fast 
hitting times, it is interesting to observe that both the quantum walk search algorithm 
on the hypercube and the glued trees graph are examples of symmetric graphs where 
the quotient graph is exponentially smaller than the original graph. In the case of the 
hypercube, the hitting time for the efi'ective walk on the quotient graph is exponentially 
smaller than the number of vertices in the original graph, and exponentially smaller than 
the classical hitting time; the same is true of the continuous-time walk on the glued-trees 
graph. It is interesting to note that in both of these cases the quotient graph is a finite 
line. This seems to suggest that if we can identify graphs which have the line as a quotient 
graph, they may be fruitful ground to look for more examples of walks with fast hitting 
times. This remains very much an open question, but it is our belief that graph symmetry 
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is of vital importance in the understanding of hitting times for quantum walks, and that 
understanding the structure of quotient graphs is the key to further progress. 
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Chapter 6 



Quantum walks in decoherence-free subspaces 



We have seen in Chapter 3, that decoherence has a detrimental effect on hitting times 
of quantum walks, especially fast hitting times. This is because the examples of fast 
hitting times that have been found in the literature are due to the symmetry of the 
graph. This would not be the case in the presence of decoherence since it would break 
this symmetry and force the walk to leave this subspace. Infinite hitting times are also 
due to symmetry and its effect of confining the walk to a subspace. We have examined 
the effect of decoherence on infinite hitting times and we have seen that it makes them 
finite by breaking the symmetry. However, if the decoherence has symmetries of its 
own, then it can have invariant subspaces and the evolution in these subspaces can be 
purely unitary. In this chapter, we explore this possibility. We derive conditions on the 
decoherence such that it has enough symmetry to preserve the evolution in some subspace. 
This can be applied to preserve subspaces that lead to fast hitting times. The concept 
of decoherence-free subspaces (DFS) has been the subject of intense research in the last 
few years. Decoherence-free subspaces are a result of the symmetry associated with the 
decoherence operators and as such are subspaces in which the evolution is purely unitary. 
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Refs. [51, 52] present overviews of the progress in this field. The conditions under which a 
given subspace becomes a DFS of the evolution is given for various different formulations 
of the evolution. Here, we first give the basic definitions and conditions of decoherence- 
free subspaces. Then we review the general theory of DFS using representation theory of 
complex algebras and apply it to quantum walks to arrive at a general condition for the 
subspace Hh to be a DFS (recall from Chapter 5, that the subspace Hh is the subspace 
of the quotient graph that gives fast hitting times for some graphs.) We also provide 
examples of decoherence which have the subspace Hh inside a decoherence-free subspace 
for both the discrete-time and continuous-time walks on the hypercube. 



6.1 Decoherence-free subspaces and quantum walks 

In order to apply the DFS formalism to the discrete-time quantum walk, consider an 
evolution given by the OSR 

V{p) = J2AipAl (6.1) 

i 

where the Kraus operators satisfy the relation Yli ^l^i = I- 

Definition 6.1. A system with a Hilbert space H is said to have a decoherence free 
subspace H' if every pure state in this subspace is invariant under the OSR i.e., V|j) G T-C 

we have 

Y,Mj)m\ = \m- (6.2) 

i 

We now state the main theorem which gives the necessary and sufficient condition for 
a subspace to be a DFS. 
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Theorem 6.2. A subspace H is a decoherence-free subspace iff the Kraus operators act 
proportional to the identity on the subspace, i.e., V \j) G H' 

Ai\j)=ci\j). (6.3) 

The above theorem states that for H' to be a DFS of the given OSR, the entire 
subspace H' must he inside a single degenerate eigenspace of each of the Kraus operators. 
This can be put in the language of the representation theory of the algebras, generated 
by the operators Ai. We state the important results in this direction. First, it can be 
shown that the Kraus operators generate a complex, associative and f-closed algebra A 
(f-closed means that if M G ^, then G A.) We also assume that the identity lies in 
this algebra. We now have the following results. 

Theorem 6.3. If A is a complex, associative, ^-closed, algebra such that I E A, then A 
has the following decomposition: 

A^ ^Inj(^Midj,C). (6.4) 
JeJ 

This states that A reduces to a direct sum of irreducible complex matrix algebras 
M{dj,C) of dimension dj with a multiplicity nj. The Hilbert space on which they 
operate can be decomposed similarly as 



(6.5) 
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Prom the theory of decoherence-free subspaces [51], the following result gives the condition 
for a subspace to be a DFS. 



Theorem 6.4. The necessary and sufficient condition for a subspace Ti! to he a decoherence- 
free subspace of the dynamics is that it must lie in the degeneracy of a single one- 
dimensional irreducible representation. 

This theorem states that in order to be a DFS, H' must lie in the part of the 
decomposition for which d^- = 1 i.e., using the decomposition in Eq. (6.5), H' C Huk- 
When 7^ 1, we obtain the more general decoherence-free subsystems. To apply this 
to quantum walks, we are interested only in decoherence-free subspaces and so we only 
need one-dimensional irreducible subalgebras. Finally, to perform computations which 
preserve the DFS, the unitary operations must commute with the algebra inside the DFS. 
This means that the general form for the operations that preserve this DFS is 



where d is the total dimension of the Hilbert space. 

We now apply this to quantum walks with decoherence and find conditions under 
which the decoherence preserves the subspace of interest-7-^jy. Firstly, using the condition 
of degeneracy in Eq. (6.3), we need that 



U G M{nK,C) ® M{d - UK, C), 



(6.6) 



(6.7) 
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where Oj is an orbit under the action of H and is a basis vector in Hh- In order to find 
a condition in terms of the algebra generated by the symmetry operators, (j{h), we first 
note that for the symmetry operators, if H' is a degenerate eigenspace i.e., 

<9)\3)=cg\j), (6.8) 

for all \ j) G H' 1 then this subspace is preserved by the unitary evolution operator U. This 
can be seen easily, for we have 

<y{9)U'\j) = U'<y{9)\j) = cgtl'lj). (6.9) 

Now, the complex algebra generated by the group of symmetry operators (denote it by 
Ah) is associative, f-closed and contains the identity. Using the above general theory, we 
can decompose this into a sum of irreducible subalgebras as 

AH=^Inj®M{dj,C). (6.10) 
JeJ' 

The condition Eq.(6.8) is equivalent to the condition that the subspace Hh niust belong 
to the degeneracy of a one-dimensional irreducible subalgebra in order to be preserved by 
the evolution and the evolution operator U must lie in the algebra that commutes with 
Ah in the DFS, i.e., if for some K' G J' for which c?'^ = 1, the subspace Hh lies in its 
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degeneracy, then U must lie in the algebra M.{n'p^,C) ® M.{d — n'j^,C). Combining the 
these conditions, we have that 

A ^ In^®C© Inj(^M{dj,C) 

JeJ-K 

Ah = In'^^C® Inj®M{dj,C) 

JeJ'-K' 

= M{nK,C)(S)M{d-nK,C), (6.11) 

where we need n'p^ = uk- Now if Hh C W^^, then the quantum walk which starts inside 
Hh will be unperturbed by decoherence. 

As an example, consider the discrete-time quantum walk on the hypercube. Consider 
a decoherence which has the following Kraus operators 

Ai = KiSii+i ® + + + (6.12) 

where z G {1, 2, . . . , n — 1}. Sij is the swap operation on the qubits numbered i and j. It 
is given by 

Sij = / (8) . . . (8)7 (8)52 (8 ®_f , (6.13) 

i— 1 n—i—1 

where 

^2 = |0)(0| ® |0)(0| + |1)(1| «) |1)(1| + |0)(1| «) |1)(0| + |1)(0| «) |0)(1|. (6.14) 
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The constants need to satisfy 

n-l 

Y^\Ki\^ = l (6.15) 

i=l 

in order to have A = I- This kind of decoherence might occur in implementations of 

quantum walks on the hypercube using n spins. The decoherence is local in the sense that 
it only affects adjacent spins. This decoherence has the subspace Hh with basis vectors 
given above as a DFS. These basis vectors and Kraus operators satisfy the condition in 
Eq. (6.7) since we have 

Ai\Oj) = Ki{Sii+i ® + 1| + |i + l){i\))\Oj) = Ki\Oj). (6.16) 

Each \Oj) is an eigenvector of eigenvalue 1 of the swap operators. 

Continuous-time quantum walks evolve by an application of a Hamiltonian as defined 
in Eq. (2.9). Decoherence in this scenario would give rise to an evolution which can be 
described using the Lindblad semigroup master equation as 

m = [H, p{t)] + J2 ai{Lip{t)Ll + L\Up{t) + p{t)Lm. (6.17) 

i 

The necessary and sufficient conditions for the existence of decoherence-free subspace for 
such an evolution is given in the following theorem [52]. 

Theorem 6.5. A subspace Ti' is a decoherence-free subspace of the above evolution if 
and only if the Lindblad operators are proportional to the identity on the subspace i.e., V 

\j) e n' 

Li\j) = di\j). (6.18) 
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An example of a decoherence model which has as its DFS the subspace T-Lh is where 
the Lindblad operators are 

Li = KiSij, (6.19) 

where Sij has the same meaning as before. Indeed, as before one can derive a sufficient 
condition on the Lindblad operators which is that they must he in the algebra generated 
by the group H. If this is true then the subspace Hh will be a DFS. 

The decoherence operators (Kraus or Lindblad) need not lie in the algebra generated 
by the symmetry operators as in the examples. The above examples consider a kind 
of decoherence that might occur in a physical implementation of a quantum walk on 
the hypercube. But in general, the decoherence operators need not lie in the algebra 
generated by the symmetry operators. The decoherence operators just need to commute 
with the symmetry operators inside the DFS. 

6.2 Discussion 

Quantum walks can have invariant subspaces if the walk has sufficient symmetries. A 
walk can be endowed with enough symmetry such that there exists a suitable invariant 
subspace, which can be used to design an algorithm. The quantum-walk based search 
algorithms exploit such invariant subspaces. If such a walk is implemented using any 
quantum system, it will have decoherence affecting it. If the decoherence has symmetries 
i.e., if the decoherence operators have degeneracies, then we can have decoherence-free 
subspaces. We have determined conditions which the decoherence operators must satisfy 
in order to have a certain subspace of interest (from an algorithmic point of view) lies 

119 



in a DFS of the dynamics. We presented examples of such a decoherence for a quantum 
walk on the hypercube. 

In general, it is difficult to know the symmetries of the decoherence operators without 
looking at a specific implementation of quantum walks. Unfortunately, there are not 
many implementations of quantum walks and the few that exist are all for the walk on a 
line. Here we present the general conditions that the decoherence must satisfy in order 
to have the invariant subspace Hh as a DFS. For a specific implementation, one can 
check if the decoherence satisfies these conditions to know if that scheme can successfully 
implement a quantum walk algorithm which exploits invariant subspaces by preserving 
these subspaces even in the presence of decoherence. 
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Chapter 7 



Conclusions 



Symmetry plays an important role in many areas of physics and especially in quantum 
mechanics. Many quantum systems have inherent symmetries which make it easier to 
analyze them. Since quantum computing involves harnessing the power of quantum 
mechanics to perform faster computations, we must take into account induced symmetries 
in order to make full use of quantum mechanics. Indeed one such example is in the theory 
of decoherence-free subspaces and subsystems. Based on the symmetries of the noise 
operators, we can determine subspaces of the Hilbert space where decoherence has no 
effect. 

In view of this, it should not be surprising if symmetry plays a very important role in 
quantum algorithms. Indeed, symmetry often contributes in an implicit way. Consider 
the Grover search algorithm where the problem is to search for a marked node among 
many unmarked nodes (for simplicity we assume that there is only one marked node) . In 
this abstract setting we can see that the labels of the unmarked nodes have no bearing 
on the problem. Hence it should have no bearing on the algorithm to search for the 
marked node. If the algorithm makes a distinction between unmarked nodes, it would 

121 



mean that we are introducing structure where it does not exist and this cannot lead to 
an optimal algorithm. Thus if we permute the labels of all the unmarked nodes in any 
way, the problem must be the same. Therefore, the symmetry group of the problem is 
Sfi-i where n is the number of nodes. Under the action of this symmetry group, we can 
easily see that there are two orbits: an equal superposition of all the unmarked nodes as 
one of them and the marked node as the other. Thus, symmetry arguments determine 
the subspace that should be the invariant subspace of the algorithm. Note that this does 
not give us an algorithm, but tells us the subspace that is important for the algorithm. 
However, in this case determining this invariant subspace is a very important step. 

One can possibly extend this to other computational problems to determine the in- 
variant subspaces relevant to them. In the case of Grover's search problem, there is no 
structure in the problem and hence every permutation of the nodes (except the marked 
node) is a symmetry. For problems with structure this will not be the case i.e., only cer- 
tain permutations of the basis elements (which need to be defined appropriately) which 
preserve the structure of the problem are symmetries. This gives a group of symmetries- 
the structure group of the problem. We can then find the orbits under the action of 
this structure group and this gives us an invariant subspace of the problem: a space 
which respects all the symmetries. Any algorithm for this problem must also respect 
these symmetries and hence must lie in this invariant subspace. Explicitly considering 
the symmetries of a problem in this way and incorporating it into the solution would lead 
to more efficient algorithms. 
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